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The use of chirally rotated boundary conditions provides a formulation of the 
£Sj Schrodinger functional that is compatible with automatic 0(a) improvement of 

Wilson fermions up to 0(a) boundary contributions. The elimination of bulk 0(a) 
effects requires the non-perturbative tuning of the critical mass and one additional 
boundary counterterm. We present the results of such a tuning in a quenched setup 
for several values of the renormalized gauge coupling, from perturbative to non- 
perturbative regimes, and for a range of lattice spacings. We also check that the 
correct boundary conditions and symmetries are restored in the continuum limit. 
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1 Introduction 



The Schrodinger functional of QCD [1 .2 ,3J is the gauge invariant functional integral for 
QCD on a hyper-cylinder where the fields satisfy periodic boundary conditions in the 
spatial directions and Dirichlet boundary conditions at the Euclidean times and T. The 
SF is non-perturbatively defined and it has been shown to be very successful when used 
as a renormalization scheme for lattice QCD. An incomplete set of results obtained using 
the SF is given by refs. 



UHfflHMIINIIIIH 



Because of the Dirichlet boundary conditions in the time direction, the SF has a spectral 
gap even at zero quark mass [2lfT3] . thus allowing to use the SF as a massless renormaliza- 
tion scheme. Moreover, due to the possibility of applying finite-size techniques, the SF is 
an ideal framework to evaluate scale-dependent quantities over a wide range of energies, 
covering the perturbative up to the non-perturbative regimes. Such a framework is needed 
when studying non-perturbative renormalization on the lattice. Another good property 
of the SF is the ability to define gauge invariant quark sources, making it possible to 
construct gauge invariant correlation functions to determine renormalization factors. 

The implementation of the SF on the lattice beyond the pure gauge theory is not a 
straightforward issue. Depending on the bulk lattice Dirac operator adopted, boundary 
terms in the lattice action may need to be added in order to recover the correct boundary 
conditions in the continuum limit. For the implementation with Wilson fermions [2113] . the 
boundary conditions arise naturally [13] and no fine-tuning is needed near the boundary. 
For Ginsparg- Wilson fermions, specific boundary terms have to be added to the lattice 
action in order to recover the correct continuum limit [T3]. A related issue is the definition 
of boundary conditions in order to achieve automatic 0(a) improvement [H] with massless 
Wilson quarks. In this case one would like to have boundary conditions in the continuum 
limit that allow Wilson fermions to maintain automatic 0(a) improvement, similarly to 
what happens with twisted mass fermions. A solution to this problem has been proposed 
by Sint [T5fTo] . called the chirally rotated Schrodinger functional (xSF). Being compatible 
with automatic 0(a) improvement, makes the %SF an ideal setup for renormalizing bare 
operators computed with Wilson twisted mass fermions at maximal twist [TTJ and it 
may have several advantages compared with the standard SF. Here we just mention the 
possibility of computing renormalization factors of four-fermion operators (like for B^) 
or of twist-2 operators (like for (x)) free of 0(a) corrections. 

The main content of this paper is the numerical investigation of the non-perturbative 
tuning of the xSF in a quenched setup. In sect. [2] we start with Wilson twisted mass 
fermions as an example of automatic 0(a) improvement. Then in sect. [3] we discuss the 
xSF in the continuum and its relationship with automatic 0(a) improvement. In sect. [4] 
we introduce the appropriate lattice action for Wilson fermions and the boundary coun- 
terterms and in sect. [5] we discuss the non-perturbative tuning of the critical mass and 
of the relevant boundary counterterm needed to obtain the correct continuum limit. In 
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sect. [6] we check that the proper symmetries and boundary conditions are recovered in 
the continuum limit. In a forthcoming paper [TH] we will study the application of the 
xSF renormalization scheme to the determination of several physically relevant quanti- 
ties. The computation of such quantities will moreover allow us to perform a continuum 
limit scaling test for cutoff effects both for small and large volume calculations. 



2 Wilson twisted mass fermions and automatic O(a) improvement 



The Wilson twisted mass (Wtm) formulation [T^PU] is a lattice action that provides a 
solid framework to perform large scale simulations with Nf = 0,2 [2T|22pi7|l23] and more 
recently Nf = 2 + 1 + 1 [21] flavors of dynamical fermions. One of its advantages is the 
automatic 0(a) improvement of physical correlation functions [H], which requires the 
non-perturbative tuning of just a single parameter: the critical mass m cr . The resulting 
lattice action is referred to as Wtm at maximal twist. Clearly it would be desirable to 
retain automatic 0(a) improvement thorugh the process of renormalizing local operators. 

We now briefly summarize how automatic 0(a) improvement works for Wilson twisted 
mass fermions in a finite volume without boundaries. The very same mechanism will be 
used in the next section to show how the %SF retains this property. The action for twisted 
mass QCD (tmQCD) in the continuum for a flavor doublet x °f fermions is 



S 



d X X(x) 7/^ + m q + ^q75T X( X ) 



(2.1) 



where m q is the so-called untwisted mass, /i q the twisted mass and r 3 is the third of the 
Pauli matrices r a . Performing the following non- anomalous change of basis 



^(x) = e^ 75r3/2 x(x), i>{x) = x(x)e iuJ ^ T3/2 , 



oj = arctan 



(2.2) 



it is easy to see that the tmQCD action (2.8) is equivalent to the standard QCD action 



for degenerate Nf = 2 fermions ip with mass M = y m q + This trivial change of basis 
becomes non-trivial once we decide to discretize the QCD action with Wilson fermions 
obtaining 

s = a4 J2x(%) v w + mo + X(x), (2-3) 



where T>w is the standard Wilson operator 

1 r 



V 



w 



(v„ + v;) 7m - «v;v M 



(2.4) 



The Wtm action (2.3) has the proper continuum limit [19] and, after tuning the bare un- 
twisted mass mo to its critical value m cr , physical correlation functions are automatically 
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0(a) improved [14^] The property of automatic 0(a) improvement of physical correla- 
tion functions is a consequence of the different transformation properties of the mass term 
and the Wilson term under vector and axial symmetries. 



We start by noting that the Wtm lattice action (2.3) is invariant under the transformation 



Til' 2 x V 1 ' 2 where 



X(x) ijsT^xix) 



is a discrete chiral transformation and 

{U (x; jj) — > W(—x — a/2; 11) 
X (x) -> ir^e^xi-x) (2.6) 
x(x)^x(-x)Hr^)e 3 ^ 2 , 

is a discrete vector transformation combined with a transformation that essentially counts 
the dimensions of the fields [H]. We consider a general multiplicatively renormalizable 
multilocal lattice field $ that is even under the transformation Til' 2 x ^ >1 ' 2 - in the following 
we will refer to even (odd) operators and correlation functions under a transformation if 
they are invariant (change sign) under that transformation. The discretization errors of 
the lattice correlation function ($) are described by the Symanzik effective theory |30j . 



The Symanzik effective action corresponding to (2.3) (with mo = m cr ) reads 

S cS = S + aSi + . . • , (2.7) 
where the target continuum theory is 

5 = Jd A x x(x) 7^ + i/iq7 5 T 3 x(x) . (2.8) 

We recall that the effective theory is constructed taking into account the symmetries of 
the lattice action [31]. This implies that the higher-dimensional correction terms in the 
effective action 

5 1 = fd 4 y C x (y) C l {y) = £ cAG/) (2-9) 

i 

are Til' 2 x T> 1,2 even. After using the equations of motion, the only operators Oi contribut- 
ing to on-shell correlation functions for vanishing untwisted quark mass are [3~Tf3"2] 

iX^^F^X, I^XX, (2-10) 

where a^ v = | [7 M ,7z,] and F^ v is the gluon field strength tensor. We observe that, even 
if both Sq and S\ are invariant under Til' 2 x D 1 ' 2 , the continuum theory is separately 



For other proofs of automatic 0(a) improvement see refs. [25 26 27 28 15 29J . 
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invariant under TZl' 2 and D 1 ' 2 , while C± is odd under both Til' 2 and V 1,2 . In the effective 
theory $ is represented by an effective field 

$ eff = $ + a$ 1 + ... , (2.11) 

where $1 represents a linear combination of 0(a) counterterms specific to the field $. A 
renormalized lattice correlation function of the field <3> to order a in the effective theory 
is then given by 

($> -> ($ )o - a Jd A y ($o£i(3/))o + a($i)o + • • • (2.12) 

where the expectation values (• • • )o are to be taken in the continuum theory with action 
S . 

If we are interested in a non-vanishing correlator in the continuum limit, ($0)0 7^ 0, $0 
must be even under both Til' 2 and V 1,2 . Because of its higher dimensionality, this implies 
that $1 is odd under V 1,2 and thus odd under Til' 2 as well. We have already noticed that 
Ci is odd under both Til' 2 and V 1,2 . We can conclude that both ($1)0 and Jd 4 y (<3>o£i(2/))o 
vanish because the continuum theory is invariant]" 3 "] under TZl' 2 . 

The key point for the absence of 0(a) terms in the Symanzik expansion of TZl' 2 even 
correlation functions is that the continuum action (2.8D is invariant under the discrete 
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chiral transformation 1Z 5 ' , while all the operators in eq. (2.10) of the Symanzik expansion 



of the lattice action are odd under the same discrete chiral 



symmetry transformation. 
Furthermore, the form of the continuum theory (2.8), with vanishing untwisted quark 
mass, that guarantees automatic 0(a) improvement is a direct consequence of the non- 
perturbative tuning of m = m cr in the lattice theory. Possible uncertainties of 0(a) in 
the determination of the critical mass m cr are proportional to XX-, hence their insertions 
in the effective theory vanish as the insertions of L\ do. 

Another way of seeing automatic 0(a) improvement is by saying that the lattice ac- 



tion (2.3) has two distinct sources of chiral symmetry breaking: the Wilson term (together 
with the critical mass m cr ) and the twisted mass term. Automatic 0(a) improvement is 
a consequence of the fact that one of the two terms (the twisted mass term in our basis) 

1 2 

retains the discrete chiral symmetry 1Z 5 ' . In the next section we show that the same idea 
applies to the %SF where now the two source of chiral symmetry breaking are the Wilson 
term and the boundary conditions satisfied by the fermion fields. 

It might come as a surprise that automatic 0(a) improvement works only for correlation 
functions that are even under 71$ 2 . To understand this we need to do a step back to our 



target continuum theory (2.8). We have shown in this section that tmQCD and QCD 



3 As usual possible contact terms in fd 4 y (3>o£i(y))o can be traded for terms with the same 
symmetry properties of ($1)0 without invalidating the proof. 



5 



are the same continuum theories written in a different fermion basis. Therefore QCD and 
tmQCD share the same symmetry properties even if the symmetry transformations take 
different forms in the different basis. In app. [A] we collect few symmetry transformations 

— 1 2 

in the twisted basis for a generic twist angle 00. The TZ 5 ' transformation applied to fermion 
fields whose target continuum theory is eq. (2.8), i.e. u = n/2, corresponds to a vector 
(flavor) transformation in the basis where the QCD action takes its standard form, i.e. 
the mass term takes its standard form. It is easy to see it considering eq. (A. 6) with 
a v = (tt, 0, 0) and 00 = or uj = ix/2. To summarize the Til' 2 transformation takes the 
form of a discrete chiral transformation, but it has the physical meaning of a vector (flavor) 

1 2 

transformation. For this reason in the following we will refer to 1Z 5 ' -even correlation 
functions as "physical" to distinguish them from the correlation functions that vanish in 
the continuum limit. 



3 Chirally rotated Schrodinger functional 



A well-known and successful non-perturbative renormalization scheme is the so-called 
Schrodinger functional (SF) scheme. The SF for QCD is the standard QCD partition func- 
tion where the fermion and gauge degrees of freedom satisfy Dirichlet boundary conditions 
at xq = and T. Periodic boundary conditions in the spatial directions are employed for 
the gauge fields, while fermion fields can be defined to be periodic up to a phase 

ip(x + Lk) = e i0k i()(x) -n<9 k <n. (3.1) 

The boundary conditions for the gauge fields are described in detail in ref. [Tf3"T] . In the 
following we will concentrate only on the fermionic fields and assume throughout that the 
boundary conditions for the gauge fields in the continuum and later on the lattice are the 
standard ones. A natural choice for Dirichlet boundary conditions (b.c.) for the fermion 
fields are the so-called standard SF b.c. [2"1fT3"] 

P + ^(x)\ XQ=0 = V P_V(x)| aio=T = fromT (3.2) 

ip(x) P_|a; 0= o = from C ip(x) P + \ Xo= t = from T and C (3.3) 

with the projectors, 

P ± = 1 (1 ± 70) , (3.4) 

and where we have specified the discrete symmetries (charge conjugation C and time 
reversal T) that relate the different boundary conditions (see app. [A] for their definitions). 
The boundary fermion fields are defined as 

P-i>(x)\ X0=0 = C(x) P + ^(x)\ X0=T = C(x) (3.5) 

$(x) P+L=o = C(x) $(x) P-\ X0=T = C'(x) . (3.6) 
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Even in the massless limit standard SF b.c. break chiral symmetry. Studying the transfor- 
mation of the SF propagator under a chiral symmetry transformation one observes that 
the SF b.c. induce a unit mass-like term at the boundaries [13] • The breaking of chiral 
symmetry by the boundary conditions implies, contrary to what happens in a finite vol- 
ume without boundaries, that Wtm at maximal twist is affected by 0(a) discretization 
errors. 

From the proof of automatic 0(a) improvement presented in sect. |2j we understand that 
the relevant property to preserve automatic 0(a) improvement is the way chiral symmetry 
is broken in the continuum theory (by a mass term or by boundary conditions) with 
respect to the way the Wilson term does it in the lattice theory. We recall that in the 

1 2 

case of Wtm without boundaries the twisted mass term is invariant under 1Z 5 ' , while the 
Wilson term is not. One might thus think to define the continuum theory with Dirichlet 
boundary conditions invariant under Til' 2 . The presence of the Wilson term at non-zero 
lattice spacing should not harm the property of automatic 0(a) improvement [15] . 

One possible solution to this problem is to mimic exactly what is done with tmQCD. 

1 2 

In the continuum, the boundary conditions preserving TZ 5 ' can be obtained from the 
homogeneous standard SF boundary conditions via the non-singlet axial transformation 



defined in eq. (2.2). If we consider a flavor doublet of fermions and we apply such a 



rotation to the quark and anti-quark fields, the boundary conditions take the form, 

Q+ X(a0 l<no=o = ° Q- x(x)\xq=t = from %/ 2 (3.7) 

X(x) Q+\ Xo =o = from C x( x ) Q-\x =t = from %/2 and C (3.8) 

with projectors 

Q± = \(l± i7o7 5 r 3 ) • (3.9) 

These are the chirally rotated b.c. pj)], which we will refer to as the %SF b.c. (see app. [A] 
for the definition of 7^/2)- The boundary fields can be defined as 

Q-X(x)L=o = C(x) Q +X (x)\ X0= t = C(x) (3.10) 

x(x)Q-\x =o = ((x) x(x)Q+\ X0 =T = ('(£)■ (3-11) 

It is important to notice that the correspondence between the SF and the %SF is analogous 
to that between QCD and tmQCD. The symmetries of the SF are the same as those of 
QCD while the symmetries of the %SF correspond to those of tmQCD at maximal twist. 
In fact, these symmetries are not different in the two formulations, they are just expressed 
in a different basis. In the continuum, the SF and the %SF have all the same symmetries. 



The transformations of eqs. (2.2) are a trivial change of basis in the continuum theory, so 
one might hope that on the lattice massless Wilson fermions with %SF b.c. will provide 
a framework for a finite volume scheme compatible with automatic 0(a) improvement. 
We have noted before though that standard SF b.c. arise naturally when performing the 
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continuum limit with massless Wilson fermions. This implies that in order to recover 
the xSF b.c. in the continuum limit additional terms have to be added to the Wilson 
action near the boundaries. This problem has been solved by Sint in [TO] using orbifolding 
techiques. Orbifolding assures that the proper b.c. are satisfied at tree-level of the lattice 
theory. This is sufficient to identify the proper terms to add to the lattice action near 
the boundaries. The study of the renormalization of the theory will teach us if additional 
terms are needed to obtain the correct continuum limit in the interacting case. In the 
next section we discuss in more details the lattice action proposed in [To] . 

We conclude this section by emphasizing that the %SF in the continuum limit is the 
same renormalization scheme as the standard SF, if the same kinematical conditions are 
chosen. It is only at non-zero lattice spacing where the two schemes differ. This property 
is important for many reasons. Here we just mention that previous results obtained with 
the standard SF can be used to check the validity of the continuum limit of the %SF. 
We will use this property in our continuum limit scaling studies in ref. [18]. Moreover, 
if one is interested in renormalizing certain operators for which the evolution with the 
renormalization scale has already been computed with the standard SF, then it is sufficient 
to compute the proper renormalization factors at the lattice spacings where the infinite 
volume operators are used. For the scale evolution the results from the standard SF can 
then be used. 



4 xSF with Wilson fermions 



The construction of a lattice action that in the continuum limit goes to QCD with %SF 
b.c. is a non-trivial task. In ref. [16J three lattice actions have been proposed which satisfy 
this property. These actions are the standard Wilson action in the bulk of the hyper- 
cylinder with three different local modifications close to the temporal boundaries. These 
modifications are necessary to obtain the %SF b.c. in the continuum limit. 

In this section we briefly discuss the lattice action that we have used in our numerical 
investigation. In what follows, we call this the Wilson %SF (W^SF) action. For additional 
details about the other two formulations see ref. [IB] . We consider a doublet of fermions 
X — (Xu,Xd) T and a lattice L 3 x [0, T] with spacing a. The W%SF action reads 

S F = a 4 Xi x ) ( p w + m ) X (x) , (4.1) 

x =0 x 

where x = (xq,x) is a point on the lattice with spatial coordinates x and temporal 
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coordinate Xq and 

(-U (x)P-x{x + a6) + (AT + i7 5 r 3 P-) X (x) if x = , 

aV wX (x) = I aD wX {x) ifO<x <T, (4.2) 

[ (K + i7 5 T 3 P+) x{x) ~ U (x - a6)tp +x (x - a6) if x = T , 

satisfying, as does the twisted mass operator, the Hermiticity property, 

r^TsPwTsr 1 ' 2 ^^. (4.3) 



_D\y is the massless Wilson operator defined in eq. (2.4) and it may be written as 

aD w xO) = -U (x)P-x(x + ad) + K X (x) - U (x - a6) j P +X {x - aQ) , (4.4) 
with K, the dimensionless time-diagonal kernel of the Wilson operator, 

K = l + jz^{lk [V* k (x) + V k (x)} -aV* k (x)V k (x)} . (4.5) 
k=i z 

The spectrum of the Hermitean lattice operator, 75r 1 ' 2 2>w, is bounded from below [16], 
as in the continuum, with a non-vanishing minimum eigenvalue which coincides with the 
one in the continuum theory in the limit a — > 0. 

We note immediately that the main difference between the W^SF and the standard SF 
is the presence of an additional term 

x(x)i-/ 5 T 3 P_ X (x)5 xofi + X (x)i>y 5 T 3 P +X (x)5 XOiT , (4.6) 

localized at the boundaries. This term is necessary but not sufficient to recover the proper 
b.c. in the continuum limit. 

To ensure the correct continuum limit, one must account for all relevant operators allowed 
by the symmetries of the action above. This means to consider operators of dimension 
four or less for the bulk action. There is one such operator, XX , and the corresponding 
counterterm is the term proportional to the critical quark mass, m cr . This is the standard 
operator that is present for all Wilson actions due to the breaking of chiral symmetry by 
the Wilson term. 

Similarly, we must include all permitted boundary operators of dimension three or less. 
Again, the one allowed operator is XX [16J, which gives rise to the following counterterm 
to the lattice action, 

5S 3 = (z f - l)a 3 J2 (xx\x =o + XX\x =t) ■ 

X 

Such an operator would be forbidden in the continuum action, but the reduced symmetries 
of the Wilson action do not allow us to exclude this operator at non-zero lattice spacing. 
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The presence of 8S3 can then be understood as necessary to restore in the continuum 
limit the symmetries broken by the Wilson term. More specifically the Wilson term, 
as the 5S 3 term, break the discrete symmetry 71$ 2 . We have noticed in sect. [2] that to 
ensure automatic 0(a) improvement we want to have a target continuum action and 
b.c. invariant under an Til' 2 transformation. The parameter Zf has to be tuned in order 
to recover the proper symmetries in the continuum limit that ensure automatic 0(a) 
improvement. Since Til' 2 is a symmetry of the massless continuum theory, which is only 
broken in the regularization procedure, Zf accounts for a finite renormalization, that is, 
it has the form, 

z f (go) = zf ) +4 ) g 2 o + 0(gt), (4.7) 
with all coefficients in the expansion being finite. However, the fact that 5Ss is not an 
irrelevant operator implies that a perturbative computation of Zf is not sufficient and that 
we then must compute the bare coupling dependence of Zf(g ) non-perturbatively. 

In perturbation theory, only the tree- level value of Zf is presently known and it takes the 
value z^ = 1. We have determined this value by a direct comparison of the free quark 
propagator in the continuum and the continuum limit of the analytical expression for 
the free lattice quark propagator [53] . The analytical expression of the lattice tree-level 



propagator for the action (4.1) is given in app. B 



Furthermore, we must also examine those irrelevant operators that in principle can lead 
to 0(a) contributions. In the bulk, there is the dimension five Sheikholeslami-Wohlert 
term, but automatic 0(a) improvement eliminates the need for this operator. Yet, there 
does remain an 0(a) contribution from the boundary due to the irrelevant dimension four 
operator [16], 

5S 4 = (d s - l)a*J2 (xjkDkX\x =o + XJkDkX\ 



|Z0 = 



where = |(V£ + V&). Such a contribution is present in all the SF formulations 
and it is not due to the particular lattice action or b.c. we have chosenp"] Given that SS4 
is an irrelevant operator, a perturbative calculation of d s is presumably sufficient and the 
expansion in powers of g\ reads 

d s (g ) = + g 2 + 0(g$ . (4.8) 



For the lattice action (4.1) d s is already needed at the tree- level of perturbation theory 
in order to remove 0(a) boundary cutoff effects. The tree- level value is for the W%SF 
action in eq. (4.1) = 1/2. We have determined this value from a numerical inspection 



of the free quark propagator on the lattice, obtained from the numerical inversion of the 
following lattice Wilson operator 

S p + 5S 3 + 5S 4 . (4.9) 
This is in complete agreement with the analytical result obtained in [T6] . 



In fact, d s plays a role that is analogous to the ct counterterm in the standard SF |31| . 
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The knowledge of guarantees boundary cutoff effects of at most 0{agl) and like in 
the standard SF, we expect that a perturbative determination of d s is enough to cancel 
the dominant 0(a) boundary effects. It goes without saying that a determination beyond 
tree- level would be very desirable. 

The important conclusion of the above discussion is that, with respect to the standard 
formulation of the SF, there is an additional boundary coefficient, Zf, which has to be de- 
termined non-perturbatively. However, this is enough to guarantee the correct continuum 
limit of the theory and bulk automatic 0(a) improvement up to boundary effects of at 
most O ( ag\ ) PH Therefore , besides the boundary improvement counterterms to the action, 
no further improvement counterterms need to be added to any 1Z 5 ' 2 -even quantity. Thus 
the action given in this section retains all the advantages of automatic 0(a) improvement. 



5 Non-perturbative tuning 

From purely theoretical considerations, we have concluded in the previous section that 
W%SF provides a suitable discretization for the %SF non-perturbative renormalization 
scheme. This is achieved in principle with the non-perturbative tuning of only two pa- 
rameters which are functions of the bare gauge coupling g^. the bare quark mass, mo, 
and the boundary coefficient, Zf. The bare quark mass needs to be tuned to its critical 
value, m cr , in order to have a massless scheme, while the tuning of the coefficient Zf to its 
critical value, zf, is required in order to recover the desired boundary conditions in the 
continuum and thus to obtain bulk automatic 0(a) improvement. 

If zf is not determined correctly, then the 7^5' 2 -symmetry would not be properly restored 
and bulk automatic 0(a) improvement would not take place. The loss of automatic 0(a) 
improvement may be the least of our worries if the continuum limit itself is compro- 
mised by incorrectly fixing zf . Clearly, a non-perturbative determination of zf is then 
mandatory. 

The non-perturbative determination of zf can be carried out by requiring suitable 7ll' 2 - 
odd correlation functions to vanish (cf. also [31]). Since these conditions are not unique, 
different determinations of zf are expected to differ by 0(a) effects, which should only 
affect 7?-5' 2 -even correlation functions up to 0(a 2 ). This is similar to what happens in 
large volume simulations with Wilson twisted mass fermions at maximal twist. The in- 
trinsic 0(a) uncertainties in the determination of the critical mass m cr only affect physical 
quantities at 0(a 2 ). 



We note that, similar to the boundary term proportional to d s , there is in the gauge action an 
improvement coefficient c± [31], which multiplies a dimension four boundary term, that is also 
only known from perturbation theory. 
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In fact both m and Zf have to be tuned non-perturbatively and simultaneously if a mass- 
less renormalization scheme with %SF boundary conditions is to be defined. In particular, 
it is very important to understand whether this 'combined' tuning is feasible at all, as 
otherwise a practical application of the %SF scheme would be rather cumbersome. 

After the proper determination of m cr , zf and the two boundary improvement coefficients 
to the action, c t and d s , automatic 0(a)-improvement is expected to hold. This means 
that without any improvement counterterm to the action in the bulk and to the fields, all 
physical quantities have leading 0(a 2 ) discretization effects. In practice, only a perturba- 
tive determination of the boundary improvement coefficients, q and d s , is available. At 
present, only the tree-level value of d s is known in perturbation theory. For q, we employ 
the 2-loop value [351136] . c t (g ) = 1 - 0.089 g$ - 0.030 gfc. 

Here we are concerned with the non-perturbative tuning of the other two coefficients, m 
and Zf. From now on, as it is usually done, all discussions will take place in terms of the 
hopping parameter k = 8+2 1 amo ■ Due to the potential complications which may arise in 
the tuning procedure, we first performed some studies at the tree-level of perturbation 
theory. We have tested several tuning strategies, and the preferred one, as it emerged from 
our tree-level investigation, was applied non-perturbatively to the interacting theory 



as will be explained in sect. 5.3 Besides the particular selection of the tuning strategy, a 



tuning condition must also be chosen. In sect. 5.2| we will describe all the tuning conditions 
we have investigated. 



5.1 Some definitions 



The non-perturbative determination of k cv and zf requires imposing conditions at non- 
zero lattice spacing that ensure the restoration of the expected symmetries in the contin- 
uum limit that are broken by the Wilson term at non-zero lattice spacing. Moreover, these 
conditions should be imposed at each lattice spacing while fixing a suitable renormalized 
quantity. In this work, we keep the renormalized SF coupling, g, fixed. This is equivalent 
to fixing the physical size of the box, L (we choose T = L). All other external parameters 
must also be held fixed. These are d s , which is set to its tree- level value, dy\ and the 
spatial momentum, p, which is set to zero. In the spatial directions, periodic boundary 
conditions up to a phase are assumed, whose phase dependence is parametrized by the 
angles 6 = (9±, 6*2, 6*3). During the tuning procedure, these angles are used in order to de- 
fine alternative tuning conditions. To be concrete, we choose the symmetric case, 9j. = 9 
(k = 1, 2, 3) and two values of 9 are used, 9 = 9^ = and 9 = 9b = 0.5. Different choices 
for 9 are used so that we can define tuning conditions that differ by O(a) The choices for 
the external parameters are summarized in tab. [TJ 

Before specifying the tuning conditions, we define the correlation functions that are needed 
for our tuning procedure. In particular, we will employ boundary to bulk correlation 
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T/L 


d s 


Pi 


P2 


P3 


0a 




x 


2/0 


1 


0.5 


0.0 


0.0 


0.0 


0.0 


0.5 


T/2 


3T/4 



Table 1 

Fixed parameters during the tuning. See the text for a discussion of these parameters. 



functions that involve the boundary at x = 0. For this purpose, we first define the 
boundary operators. A definition of the fermion fields at the boundary (x = 0) consistent 
with gauge invariance is given by 



C(x) = C/ (0, x) X (a, x) C(x) = x(a, x)£/ (0, x)* 
The boundary interpolating fields at xq = are given by 



(5.1) 



<^ = a 6 £c(y)r 5 Q ± C(z) 



(5.2) 



In this expression, Q± are the %SF projectors defined in sect, I] Tg contains the flavor 
and Dirac structure of an operator of type O. Specifically for a pseudo-scalar density and 
an axial- vector current, we have 



^ = a 6 EC(y)75 T Q ± C(* 



^ ± = « 6 EC(y)7,75^±C(z). 

y,z z 



(5.3a) 
(5.3b) 



We note that there is a little difference with respect to the SF formulation, where the 
projectors are included in the definition of the boundary fields, (, Here we insert the 
projectors directly in the correlations functions to have the freedom to consider correlation 
functions with the "wrong" projectors [34]. These correlation functions ought to vanish in 
the continuum limit if the correct %SF b.c. are recovered, and this is confirmed numerically 
(see sect. 



6.3). 



Considering the previous definitions of the boundary interpolating fields, we may introduce 
now our notation for the boundary to bulk correlation functions. Given a bulk operator, 
X a (x), the type of correlation functions that we consider here are the following, 



sSLM 



a 

L 3 
a 3 



(5.4a) 
(5.4b) 
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For the tuning we used only the particular cases, 



a 3 



Af^^^D^^)- (5.5a) 
8t(^9) = 4EW®. (5.5b) 

X 



X 



We denote the correlation functions in this work with gx ; i n order to distinguish them 
from the corresponding correlation functions in the standard SF usually denoted with fx- 
One may interprete the gx as correlators in the x-basis, while the fx refer to the standard 
basis. The superscripts, a, b, denote the flavor index. The subscripts X, in gx, indicate 
the corresponding operator inserted in the bulk of the lattice. Three bulk operators are 
considered here: the pseudo-scalar density, P a (x), and the axial- vector and vector cur- 
rents, A a (x), VHx). Depending on what %SF projector, Q±, is chosen in the correlation 
functions, we have the corresponding subscript ±. Due to the particular %SF boundary 



conditions, cf. Eq. (3.7)-(3.8), all correlation functions defined through Q + at x$ = 
should vanish in the continuum limit. The same holds for Q_ at xq = T. However, we do 
not consider here such correlation functions. Therefore, such kind of correlation functions 
will be used only later on to perform checks on the recovery of the correct b.c. in the 
continuum limit. For the tuning conditions we consider only correlation functions defined 
through Q_. 

As a last consideration before going into the details of the particular tuning conditions, 
we also need to define the correlation function, 

The notation is the following. Let us consider the improved axial current, 

(A I )Xx)=Aftx) + ac A d lt P*{x), (5.7) 
where the derivative on the lattice, d^, is defined to be the symmetric derivative, 

K = \ (9*, + 9,) , (5.8) 

with the standard definition of the partial derivatives on the lattice. The correlation 
function (gi) A ± (xo, 0) is defined as 

(<7i)t ± (*o,0) = -^E(W^I) • (5-9) 
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This is just the equivalent of Eq. (5.5b ), where the expression of the improved axial current 
is used instead of the unimproved one. Eq. (5.9) may be rewritten in terms of Eq. (5.5a) 
and Eq. (5.5b) as follows, 



(9i)t + Oo, 9) = gT + (x , 9) + c A s(x„, 9) 



Jib 



s(x ,9) = ad fM gf ± (x ,9). (5.10) 



By substitution of Eq. (5.10) into Eq. (5.6), the last can be cast in a more explicit manner, 

Gab 
i 



r A ,(x ,y 0] e,6')=gl b ,(x ,e) 



d„gt{x ,9) 



(5.11) 



d»gt(y ,9)-d IM gpAy ,9>) 



This expression is independent of the improvement coefficient of the axial current, ca, and 
it indicates that, G^ ± (xo, yo; 9, 0') equals g A ^ ± (xo, 9) up to cutoff effects of leading 0(a). 

The correlation function G^ ± (xo, yo; 9, 9') is useful because it allows us to define yet an- 
other tuning condition differing by 0(a) effects. In particular, in the way it is constructed, 
this correlation function would not depend on the cutoff effects that in the standard setup 
are removed by employing a non-perturbatively tuned ca improvement coefficient. 



5.2 Tuning conditions 



As already anticipated, imposing distinct symmetry restoration conditions would give rise 
to different values of K cr and zf due to cutoff effects. Therefore, it is important to study 
the sensitivity of k and Zf to the particular definitions in order to better understand 
the intrinsic uncertainty in the determination of these counterterms. Our exploratory 
studies [37] have shown that the tuning of n does not pose any special problem in the 
xSF setup with respect to the standard formulation with Wilson fermions, thus we have 
decided to investigate only one tuning condition for k. On the contrary, as it will be 
shown later, there is a rather big sensitivity on the choice of the tuning condition used to 
determine zj . Therefore, we have concentrated our efforts in the investigation of different 
tuning conditions for Zf, where we have studied seven different possibilities, which we 
number from (1) to (7). Note that even if different tuning conditions have been used to 
define zf, the tuning strategy we have employed is the same for all the tuning conditions 



(cf. sect. 5.3) 



To tune k to its critical value we adopt the standard procedure of imposing a vanishing 
PCAC mass. To be concrete, it is defined here as, 



dog a _( x o, 9a) 
^gjjxoj, 



m PCAC = iTTT7= n \ ■ ( 5 - 12 ) 
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To tune Zf we require a 7^5 ,2 -odd correlation function to vanish. The correlation functions 



(from (1) to (7)) that we have used to tune Zf are the following, 



(1) 


= flA _ 


= 9a _(x ,9 a ) , 


( 2 ) 


/ 

= #A - 


„11 (-= n \ 


(3) 


= n diS 


= fl'Ao- — fl'Ao- ' 


( 4 J 




1 3 

— 7T 12 (~ a \ 
= o <L 9v k _{Xo,VA) 






6 k=l 


(5) 


= 9v k _ 








6 k=i 


(6) 


- 9v k _ 


= 9v k _ - g'v k _ , 


(7) 


= G A 


= G Ao _(x ,y ; 9a, Ob 



(5.13a) 
(5.13b) 
(5.13c) 

(5.13d) 

(5.13e) 

(5.13f) 
(5.13g) 



The values of the parameters used in the eq. (5.12) and eq. (5.13) are given in tab. [TJ In 
all these conditions, the particular combinations of interpolating fields with their corre- 
sponding Dirac and flavor indices are chosen such that the resulting correlation function 



is 7^5 ,2 -odd and non-vanishing by definition at non-zero lattice spacing. That is, the cor- 
relation functions do not violate any symmetry of the lattice theory. Eq. (5.12) and the 
conditions (l)-(6) are obtained directly from the definitions in Eq. (5.5) with the cor- 
responding substitutions. The condition (7) is obtained from Eq. (5.6), also with the 
corresponding substitutions. 



All the different tuning conditions introduced in this section allow us to study the de- 
pendence of the 'physical' correlation functions on the different values of zf obtained. 
In addition, having a number of tuning conditions at our disposal enables us to test the 
universality of the continuum limit. 

In this work, as is the usual choice in SF schemes, we have defined the correlation functions 
in the middle of the time-extent of the lattice, Xq = T/2. The only exception is the 
condition (7), which involves two time slices. There the choice is Xq = T/2 and yo = 3T/4. 
The reason for these choices is to stay as far away as possible from the boundaries, thus 
avoiding boundary effects. In case of (7), the condition was to maximize the distance 
between the two time-slices, while still staying away from the boundaries. 

Before turning to our tuning strategy, we need to make one last remark regarding our 
particular choices of tuning conditions. In order to restore the symmetries of the theory, 
we impose that the different symmetry-violating correlation functions vanish at non-zero 
lattice spacing. However, to remove cutoff effects at tree-level, a better choice would be to 
force the corresponding correlation function to take its tree-level value at non-zero lattice 
spacing. From our studies at tree-level [33j , we have seen that such effects are very small, 
well below our statistical accuracy, and do not change our final results. 
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5.3 Tuning strategy 



To check the viability of the tuning strategy for k and Zf non-perturbatively, we have 
performed a tuning at three values of the renormalization scale \i = 1/L using only 
method (1). This situation corresponds to the results presented in [37J, where we first 
explained our tuning strategy at the non-perturbative level. This separate analysis has 
been useful to check the tuning procedure. Results obtained by this analysis are labelled 
here as obtained with method (1*). The three scales correspond to a hadronic (g 2 fixed 
with L = 1.436 ro), an intermediate (g 2 = 2.4484) and a perturbative (g 2 = 0.9944) scale. 



In sect. 5.4 we will present results obtained following the strategy presented in this section 



for all the other tuning conditions defined in sect. 5.2 



The same tuning strategy has been used for each value of ft and the corresponding value 



of L/a. The values of (5 used are given in tabs. C.l C.3 and are taken from [5j. The tuning 
is performed in several steps. 

• We calculate ampcAc and gx _ at four values of Zf, and for each value of Zf, we use 
four values of k, thus giving 16 pairs of n and Zf. This allows us to determine g& _ as 
a function of am PC Ac f° r eac h value of Zf, as illustrated in fig. [TJ 

• For each value of Zf, we perform a linear interpolation of g^ _ in terms of ampcAc to 
the point ampcAc = 0. This determines the values of gp^ _ at ampcAc = 0, denoted 
9%o-i f° r eacn °f the four values of Zf, as shown in fig. [I] as the filled symbols. 

• We now interpolate these values of g\ o _ as a function of Zf to the point of vanishing 
#A -> ^ nus giving us the critical value zj, as shown in fig. [2| 



All the numerical data for these intermediate steps can be found in ref. [33]. Next we 
determine k ct . 

• Using the same 16 pairs of k and Zf, we calculate ampcAc as a function of k for each 
Zf. This is shown in fig. [3j Note that ampcAc has a very mild dependence on Zf, so the 
four curves at fixed Zf are nearly indistinguishable. Interpolating ampcAc in K to the 
point of vanishing PCAC mass, k*, we obtain the values of k* at each Zf. The resulting 
values of function of Zf are shown in fig. Hi 

• We now interpolate these results in Zf to the previously determined value of zj v , thus 
determining the value of k ct . 



All the tuning results can be read off from tab. C.6 and tab. C.7 



A key observation of this work is the mild dependence of ampcAc on Zf, at least in the 
region near k ci and zj . This can be easily seen in fig. |3j The consequence of this is clear 
in fig. |4j the determination of k ct has a weak dependence on zf and the errors of both of 
them are relatively independent. 
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Fig. 1. 5a - vs - am PCAC at four values of zj (open symbols). All fits are linear in ampcAC- The 
values of <?a - & t a^-PCAC = 0, denoted g% _, are also plotted (filled symbols). 
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Fig. 2. <7 Ao _ vs. Zf (red). The fit is linear in Zf. The value of z?, z f{9*A _ = 0)> is also plotted 
(blue). 

5.^ Tuning results 



With the strategy discussed in sect. 5.3, we have performed the tuning of k and Zf using 



t^pcac m Eq. (5.12) and the conditions (l)-(7) defined in Eq. (5.13). The tuning has been 



performed for five fixed values of the renormalized gauge coupling, g(L), which correspond 
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Fig. 3. ampcAC vs - K a t four values of zj (open symbols). All fits are linear in k. The values of 
k at ampcAC = 0> denoted k* , are also plotted (filled symbols). 
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Fig. 4. k* vs. Zf (red). The fit is linear in Zf. The value of z™ (cf. fig. [2]) and of n cr , defined as 
k*(zj = Zf T ), are also plotted (blue). 

to five values of the physical energy scale, 1/L. In particular, the physical scale ranges 
from the purely non-perturbative to the perturbative regime and 14 values of (3 have 
been considered within that range. For better clarity, the tuning points are summarized 
in tab. [2] The notation in this table is the following. 'Scale' refers to the physical scale, 
namely, the fixed value of the renormalized gauge coupling. We have denoted the five 
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16 


12.0000 


2P 




16 


24.0000 


PP 



Table 2 

Summary of all the points where the tuning was performed. 

scales as 'NP', T, 'P', '2P' and 'PP', from the hadronic to the most perturbative scale. 
NP corresponds to L = 1.436 r , I to g 2 = 2.4484 and P to g 2 = 0.9944. For 2P and 
PP we have not determined the gauge coupling explicitly. These two scales have been 
considered in order to study the dependence of Zf on g Q , for small values of go, thinking 
of a future perturbative determination of zj, for which the knowledge of the renormalized 
coupling is not necessary. The results obtained for n cr and zf using the strategy outlined 



in sect. 5.3 for all the tuning methods are summarized in tab. C.6 and tab. C.7 for k ct 
and zf ', respectively. We also present tables showing the values used for Zf and k, at each 
of the points where the tuning was performed. In these same tables, the column labelled 
'-Wconf' represents the number of configurations used in the computation of all observables 



at the corresponding point. These are tab. C.l to tab. C.5 



During the tuning, we have used several combinations of k and Zf. As indicated in sect. 5.3 
the usual choice is to use 4 values of k and 4 values of Zf. However, there are cases where 
we have used, instead, 5 values of Zf and/or 2 values of k. In particular, we have used 2 
values of k at all the (3 values where we also performed the separate tuning using (1*). 
The reason is that, relying on the very weak dependence of k and Zf on each other, as 
it was shown in figs. [3] to |4j we expected the value of n CT not to change appreciably even 
if zj would vary visibly from one method to another. Therefore, we considered that the 
value of /t cr obtained using (1*) was already a very accurate guess on where the critical 
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value of k should be using all the other conditions. In fact, these expectations were later 
confirmed from our results of k ct , which, from one method to another, are the same within 
statistical errors. Actually, in most cases K cr did not change in any digit between any of 
the methods employed in the determination of zj (cf. tab. C.6). 



On the contrary, changes in zj among the different methods are particularly manifest 
This may be seen better in tab 



C.7 



zy does not 



Here we can see how, in most cases 
agree within errors from one method to another. This behavior becomes stronger at lower 
energies, i.e. for decreasing values of 0. Even if zj does not agree from one method to 
another, the differences are expected to be only 0(a) discretization effects and, as such, 
should vanish in the continuum limit linearly in the lattice spacing. 



In order to check this expectation, we have performed the continuum limit of differences 
in zj, as determined from different methods, at the lowest value of the renormalization 
scale (cf. also [21]). In particular, the data correspond to the differences 



Azf(m) 



(1) 



Zj [m) 



m 



2,3,4,5,6,7. 



(5.14) 



The data for Azf{m) are presented in tab. C.8 for all the methods and the corresponding 



fits, linear in a/L, are given in tab. C.9 where the point L/a = 8 has been excluded 
from all the fits. The data for Az c j(m), together with the extrapolation to the continuum 
limit are plotted in fig. |5j From this analysis we can conclude that the differences in zj 
from different methods are only cutoff effects of 0(a), as expected, which vanish in the 
continuum limit. This result may be considered as an additional test of the universality 
of the continuum limit. Moreover, discrepancies of 0(a) between different values of zf 
should affect physical observables at 0(a 2 ) at most. This expectation will be confirmed 
in the following section, where we analyze the dependence of several quantities on the 
particular tuning condition. 



5.5 Conclusions on the tuning 



We have presented the results of the non-perturbative tuning of k and Zf for the xSF at 
several physical scales, for a range of lattice spacings, and using 7 different definitions of 
zj . Our results demonstrate that the tuning of these two coefficients is indeed feasible, 
at least in the quenched approximation. Moreover, we observe that the tuning of Zf 
and k are nearly independent. This observation is important, having in mind dynamical 
fermion simulations; if this behaviour persists with dynamical calculations, it may ease the 
numerical effort necessary to perform the tuning, thus reducing the number of required 
simulations. We have also shown that even if zj differs from one method to another at 
finite values of the cutoff, such differences are only 0(a) discretization effects, as expected 
theoretically. These discrepancies vanish in the continuum limit, which itself provides 
numerical evidence of the universality of the continuum limit. 
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minus z 



as determined from different methods (cf. eq. (5.14)). The 
y r (m), determined from any other method m = 2 . . . 7. The 



data are presented in tab. C.8 All extrapolations to the continuum limit are linear in a/L and 
the point L/a = 8 is excluded from all the fits. The results from the fits are presented in tab. C.9 
The points from the different methods have been plotted slightly displaced from each other. 

6 Scaling studies and universality of the continuum limit 



In the present section, we show the results of our scaling analysis of several correlation 
functions that have been computed using the values of the critical parameters, k ct and 
zf, as determined from each of the 7 tuning conditions defined in sect. |5.2| We have 



carried out these studies at all the (3 values at which the tuning has been performed. In 
sect. 5.4 we have shown that different definitions of zf lead to critical values of Zf that 
differ from each other by cutoff effects of O(a). With the scaling study presented here, we 
demonstrate that these discrepancies in zf do not influence the continuum limit value of 
physical observables. This is a very important result since the continuum limit must be 
independent of the particular definition of the critical parameters. Furthermore, we will 
show that physically relevant quantities, when determined from the different values of 
zj T , agree within statistical errors already at non-zero lattice spacing, even at the coarsest 
lattices. This agreement holds even at the matching scale with the hadronic scheme, where 
cutoff effects are expected to be largest. Indeed, the agreement at non-zero lattice spacing 
indicates that the discretization effects induced by the O(a) uncertainties in zf are very 
small. 

In order to analyze the different co rrela tion functions, we have classified them in three 



types. The first, discussed in sect. |6.1| are those 7?.5' 2 -even correlation functions that 



have a non- vanishing continuum limit, which we refer to as 'physical'. These are the only 
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quantities which are expected to be automatic 0(a)-improved (up to boundary effects), 
provi ded k and Zf are correctly tuned to their critical values. The second kind, detailed in 



6.2 



are those 7?.5' 2 -even correlation functions which vanish in the continuum limit, 



sect. 

if the correct %SF b.c. (3.7) are recovered in the continuum limit. The last type, sect. 6.3 
are 7^5 ,2 -odd quantities. They should vanish in the continuum limit up to 0(a) cutoff 
effects if 7^5 ,2 -symmetry is restored in the continuum limit. 



These 3 kinds of correlation functions have been obtained from the definitions of the 



boundary to bulk correlation functions given in eq. (5.5). For unexplained notations in 
this section, the reader is referred to sect. 15.11 With suitable combinations of the Dirac 



odd under TZl' 2 transformations. All correlation functions labelled with + should vanish 
in the continuum limit if the correct boundary conditions are recovered, independently of 
symmetry considerations. On the contrary, all those labelled with — are a priori different 
from zero, unless some symmetry requires them to vanish. 



and flavor indices in eq. (5.5), we can define correlation functions which are either even or 



Boundary to bulk correlation functions are normalized in a standard fashion with certain 
boundary to boundary correlation functions in order to cancel the renormalization of the 
boundary quark fields. In particular in this work, only one such boundary to boundary 
correlation function is considered. It is the equivalent of fi [38] in the standard SF and it 
is defined as, 

gf(e) = ~(V' a + V b _). (6.1) 



Note that the combination of signs in Eq. ( 6.1[ ) is the only possibility for gf b not to vanish 
in the continuum limit, according to the boundary conditions satisfied in the continuum. 



6.1 TZ 5 ' -even correlation functions 



1 2 

For 14 values of /3 and several kinematic conditions, we have analyzed the 1Z 5 ' -even corre- 
lation functions gp\, gtf o _ and gl 1 as defined in eq. (5.5) and in eq. (6.1). The correlation 
functions determined at the values of n CT and obtained from the tuning conditions (1) 
to (7) for all the /3 values and all the renormalization scales considered are collected in 
tabs.[Cl0HCl9j In tabs.[Clo](ai4]we give results at 6 = (0, 0, 0), while in tabs. |C35]|Cl9 



we collect results for 6 = (0.5,0.5,0.5). 



The determination of all these quantities from the methods (1) to (7) has been performed 
via interpolations to n cr and zf. To check that the interpolation does not introduce 
additional systematic errors, we have computed observables directly at k ct and zf, without 
performing any interpolation from the first tuning method (cf. method (1*) in tabs. C.6 
and C.7). These data, computed at 6 = (0.5,0.5,0.5), are presented in tabs. C.15 C.17 



We have checked for this particular tuning method that the quantities obtained via an 
interpolation of the data do agree within errors with those obtained by means of new 
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computations performed directly at the critical values of k and Zf, denoted as method 
(1*). We therefore believe that the interpolation that we perform here does not induce 
additional errors. 

From this analysis we have found that for all the 14 values of j3, all 7^5 ,2 -even quantities 
with a non-vanishing continuum limit do not depend on the definition of zf within statis- 
tical errors. This holds for any of the values of the kinematical parameters that we have 
investigated. In fig. |6]we show, as an example, the dependence on the tuning conditions for 
gp 1 at the matching scale, L = 1.436 Tq, and for 6 = (0,0,0). Similar results are obtained 
for the other correlation functions. 

The independence at each value of the lattice spacing on the tuning conditions, i.e. on 
the particular definition of zf adopted is reassuring: no large cutoff effects are introduced 
depending on the choice of the tuning condition. The continuum limit of renormalized 
quantities and their dependence on the tuning conditions will be discussed in the com- 
panion paper [T8] . 



6.2 Recovery of the x^F boundary conditions 



We present here results for the 1Z 5 ' -even correlation functions which should vanish in the 



continuum limit, due to the particular form of the xSF boundary conditions in eqs. (3.7 



3.8) satisfied by the fermion fields in the continuum limit. We collect our results for 
the gp l + and g-^ 0+ correlation functions at all the /3 values and all the renormalization 



conditions in tabs. CIO C.19 
for d = (0.5,0.5,0.5) in tabs. 



We provide results for 6 = (0, 0, 0) in tabs. C.10 C.14 and 



C.15 C.19 



We found that for all the tuning conditions gp 1 and g^ o+ vanish in the continuum limit, 
signalling that the proper xSF b.c. are recovered. We have observed that the values of gp 1 
and gy o+ would change substantially if we vary the values of zf within their statistical 
errors. If we want to perform the continuum limit at fixed "physical" scale we have to 
take this variation into accountj^J We do this by propagating the statistical error of zf 
when interpolating gp 1 ^ and g^ 0+ in Zf. We note that gp Y + and g^ 0+ do not show the same 
behaviour for variations of K cr . 

We show in fig. [7] the behaviour of gp x + (0 = 0) towards the continuum limit. This is 
an example of a quantity that should vanish in the continuum limit if the proper b.c. 
conditions are recovered. Similar plots can be obtained for other quantities looking at 
the data in tabs. |C.10||C.19| where the first error is statistical while the second, where 



available, contains the propagation of the error of zf. As discussed in the previous and 
in the next section, we stress that this is the only case where we have observed that the 



We acknowledge a very important discussion with S.Sint on this point. 
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Fig. 6. Comparison between different tuning conditions using the 1Z 5 ' -even quantity gp_ 



The 

scale is NP and 6 = (0,0,0). Data for all methods (1) to (7) are presented (cf. tab. C.10). No 
continuum limit is performed, since this quantity takes a finite value in the continuum limit 
only after renormalization. The purpose of the plot is to compare the results from the different 
tuning conditions at non-zero lattice spacing and a full agreement can be seen. The data from 
the different methods have been plotted slightly displaced from each other. 

statistical error of zf shows significant effects in the correlation functions. 



6.3 Til' 2 '-odd correlation functions 



Amongst all the Tll ,2 -odd correlation functions, we consider correlation functions that van- 
ish only because of symmetry considerations, i.e. we do not consider correlation functions 
vanishing in the continuum limit because of the recovery of the %SF boundary conditions. 
Examples of such correlation functions are gj^_ and (Jy. ■ It is important to note that the 
same correlation functions have been used in the determination of zj . When studying the 
dependence on the tuning condition towards the continuum limit of a given correlation 
function we have obviously excluded the value of zf obtained imposing that the same 
correlation function vanishes at each value of the lattice spacing. We have though studied 
the dependence of all the other tuning conditions towards the continuum limit. 
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Fig. 7. Approach to the continuum limit of the 1Z 5 ' -even quantity gp . The scale is NP and 



= (0,0,0). Data for methods (1) to (6) are presented (cf. tab. C.10). g^ + is plotted here as a 
function of a/L. The error bars contain the statistical errors of the correlation functions together 
with the propagation of the statistical error on . 



The numerical results obtained at 8 = (0,0,0) are presented in tabs. C.20 C.24, and the 
results for 6 = (0.5, 0.5, 0.5) are given in tabs. |C25HC29 



Although results from different definitions of zf do not coincide at non-zero lattice spac- 

1 2 

ing, all the TZ 5 ' -odd correlation functions vanish in the continuum limit independently of 

m 1 2 

the tuning condition adopted. This is strong evidence that 1Z 5 ' -symmetry is restored in 
the continuum limit. 

As an example we show in the plot of fig. [8] the continuum limit approach for gj^ o _ with 
6 = (0,0,0) and in the plot of fig. [9] the continuum limit of g~y k _ for 6 = (0.5,0.5,0.5). 
These plots correspond to our most non-perturbative point, L = 1.436 r , which is the 
case where the cutoff effects are expected to be stronger. The data show a linear behaviour 
in a/L and they have been fitted with a linear fit in a/L, 



f = b + h 



(6.2) 



We have not consi dered the point L/a = 8 in th e fits. The results from the fits are 
summarized in tab. C.30 for q\\ and in tab. C.31 



for g\j . As anticipated all the 7Z 5 



1.2 
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Fig. 8. Extrapolation to the continuum limit of the T^^-odd quantity g^ . The scale is NP and 
= (0,0,0). Data for methods (2) to (7) are presented (cf. tab. C.20). The fits are all linear 
in a/L (cf. tab. C.30). The point L/a = 8 has been excluded from all the fits. We show the 



data for all tuning conditions except condition (1) because it corresponds to imposing g^ _ = 0- 
The (vanishing) continuum limit values obtained from the different tuning conditions have been 
plotted slightly displaced from each other. 

odd correlation functions vanish within errors in the continuum limit. Not being automatic 
0(a) improved they scale linearly in a/L independently of the tuning condition adopted. 



6.4 Conclusions on the scaling analysis 



From the analysis performed at several (3 values at fixed renormalization scale, we learned 
that all quantities that should vanish in the continuum limit, either by boundary condi- 
tions or symmetry restoration, have a visible dependence on the tuning condition and the 
corresponding values of zf. 

The 7?-5 ,2 -even correlation functions that should vanish in the continuum limit, if the 
correct boundary conditions are recovered, show the strongest dependence on the way zj 
has been determined. In particular we have observed that to properly study the continuum 
limit we cannot simply take the central vale of zj, but we need to propagate the statistical 
error of z c I into the error for the correlation functions, i.e. the statistical fluctuations in 
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0.02 0.04 0.06 0.08 



0.1 0.12 



Fig. 9. Extrapolation to the continuum limit of the 7^5' 2 -odd quantity . The scale is NP 
and = (0.5,0.5,0.5). Data for methods (1) to (4) and (6) to (7) are presented (cf. tab. C.25). 
The fits are all linear in a/L (cf. tab. C.31). The point L/a = 8 has been excluded from all the 
fits. We show the data for all tuning conditions except condition (5) because it corresponds to 
imposing g~y = 0. 



zf are visible in the final values for the correlation functions. If one takes the statistical 
error of zf into account in the determination of the error of the correlation functions, we 
observe, for all the tuning conditions, that the proper boundary conditions are recovered 
in the continuum limit within errors. Additionally this result is confirmed for all the tuning 
conditions adopted. 

m 1 2 

This observation is in stark contrast with the 1Z 5 ' -even correlation functions containing 
the non-vanishing component of the boundary fermion fields. In this case not only the 
statistical error on zf is irrelevant, but even the choice of the tuning condition for all 
practical purposes does not change the final values for the correlation functions even 
at non-zero lattice spacings. Although we can not make a general statement about all 
possible physical observables, our expectation is that other physical quantities, different 
from the ones studied in this work, will behave in the same manner. This is reassuring for 
further uses of the %SF for physical applications with dynamical fermions. 

1 2 

In the case of the 1Z 5 ' -odd orrelation functions, we have shown that they vanish in the 
continuum limit with leading O(a) discretization effects for all the tuning conditions, as 
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expected. This is strong numerical evidence of the restoration of 1Z 5 ' -symmetry in the 
continuum limit and the independence of the continuum limit on the particular tuning 
condition. This result is also an indication that the correct boundary conditions are re- 
covered in the continuum limit. 



7 Conclusions 



Large scale simulations with Wilson twisted mass fermions at maximal twist need a renor- 
malization scheme that preserves the property of automatic 0(a) improvement. The RI- 
MOM scheme is consistent with this requirement and results for Nf = 2 dynamical 
fermions have been obtained in ref. [39] and preliminary Nf = 4 results have been pre- 
sented in ref. [10] • Recently the x-space scheme has been tested on the Nf = 2 twisted 
mass ensembles [41] . One of the problems with these schemes is that it is difficult to 
cover the large range of scales necessary to bridge the perturbative and non-perturbative 
regimes. 

Finite volume schemes have been developed to tackle this problem. These schemes can be 
used to perform the continuum limit of step-scaling functions and to carry out the non- 
perturbative renormalization, especially for scale-dependent quantities. Thus it is very 
desirable to have the possibility to use finite volume renormalization schemes together 
with Wilson twisted mass fermions. In this work we have made a detailed study of the 
xSF scheme proposed in ref. [IE] with quenched Wilson fermions. Bulk automatic 0(a) 
improvement is achieved with a single non-perturbative tuning of the parameter Zf to its 
critical value Zf T , together with the usual tuning of n CT . We have proposed a tuning strategy 
for the simultaneous determination of zj and k cv and we have perfomed a feasibility 
study of several tuning conditions, showing that the tuning is affordable with the current 
computer resources and it does not pose any specific problem. 

The study presented in this work has shown that, as expected, unphysical correlation 
functions vanish in the continuum limit with 0(a) corrections independently of the tun- 
ing condition adopted and that the proper boundary conditions are recovered in the 
continuum limit. In addition physical correlation functions are rather insensitive to the 
method employed to determine n cr and zj . This is a very promising result since it will 
ease the computational effort required to go beyond the quenched approximation. 

In a companion paper [12] we will further investigate the continuum limit scaling behaviour 
of physical quantities that are renormalized through the xSF scheme employing the results 
of K cr and zf determined in this work. 
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A Symmetries in the twisted basis 



We list in this appendix the form of a few symmetry tranfosrmations in the twisted basis. 
The twisted basis {x , X.} an d the standard basis {ip , ip} are related by the non-anomalous 
axial transformation 

ijj(x) = exp (i | 75 r 3 ^j x(x) , ip{x) = x(x) exp (i | 75 r 3 \ . (A.l) 

Charge conjugation C is invariant under this basis transformation 

(U(x;n) -> U(x;n)*, 
X (x) C- l x{x) T , (A.2) 
X(x) -> -x(x) T C, 

where C satisfies 

-y^ = C^C- 1 , 75 = C 75 C- 1 . (A.3) 

Parity and time reversal are affected by the rotation and take a different form in the 
twisted basis. In this basis they are denoted and 7I>, respectively, and have the expres- 
sions 

!U {x ,x) U (x ,-x), U k (x ,x) -»■ Ul(x ,-x-k) 
X{x ,x) -> 70 exp(iw7 5 r 3 ) x(^o, -x) (A.4) 
X(x ,x) -> x{xo,-x) exp (2 w 75 r 3 ) 70 



7; 



-x 



?) , U k (x ,x) -> U k (-x ,x) 



x(xo,x) -> 27075 exp (2 w 75 r 3 ) x(-zo,£) 



(A.5) 



X(x ,x) -> x(-x ,x) exp (2 w 75 r 3 ) 2 70 75- 

To obtain the form of a parity or a time-reversal transformation in the standard basis, V 
or T, it is sufficient to set u = in eqs (A.4) and (A.5). 

In the twisted basis the SU(2) vector transformation takes the form 



SU(2) U : 



X(x) -> exp 

x( x ) -> 



2 75^ 



cxi> (^ i" a ) exp (z I 75 r 3 ) xfa), 



exp | 75 r 3 ) exp f — i ^ r a ) exp 



(A.6) 
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B The free lattice quark propagator for the %SF with Wilson fermions 



In this appendix we obtain the analytical expression of the quark propagator S (x, y) on 



the lattice at tree-level of perturbation theory for the action in eq. (4.1). The derivation 



of the propagator is rather standard (see for example ref. [12]) once we have the exact 
b.c. satisfied by the fermion fields at finite lattice spacing a. They can be obtained as a 
spinoff of the orbifold construnctions [16] and read 



(B.l) 



Q + (l - ^d*) S(x,y)\ xo=0 = Q- (l + ^o] Si,,, y)|„,^. 0. 
The problem we want to solve is then 

(Dyr + mo) S (x,y) = a' 4 



(B-2) 



with b.c. (B.l) where 6 xy is the dimensionless Kronecker delta and D w denotes the 



massless Wilson operator (4.4). 



The result of the calculation can be cast in the form 

S (x, y) = [Dy? + m ) G (x, y) 



with 



(B.3) 
(B.4) 



where 



G(x ,y ;p) 



1 



2u(p+)A(p+)D(pi 



e -w(p+)|x -yo| _ e -u(p^)(2(T+a)-\xo-yo\) 



- i^ 5 r 3 e-^ ){xo+yo+a) 

_|_ ^ 7o75T 3 e -w(p + )(2(T+a)-(x +j/o+a)) 



(B.5) 



Here we list all the definitions of the functions useful for the determination of the propa- 
gator. 



pt=p fl ±9p/L, #o = 0. 



= - sin (apf) , M(p ± ) = m + ]-ap 
a — 



±2 



sin 



ap 



(B.6) 
(B.7) 



The function u>(p*~), such that po = iu>(p*~), is given by 



sinh 



if) 



a*i% 2 + {A{F)-\ 



2^1/2 



(B.8) 
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and 



(^{p^) = — ipQ = - sinh 
a 



au(p i 



D (p±) = 1 + e ~2^)(T+a) 



(B.9) 
(B.10) 
(B.H) 



The analytical expression of the lattice quark propagator (B.3) has been numerically 



cross-checked with the propagator obtained from a numerical inversion of the free lattice 



Dirac operator given in Eq. (4.2) and also with the corresponding propagator in 
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C Tables of numerical results 

C.l Tuning data 



Guess values for the tuning 
Hadronic scale: L = 1.436 ro 


L/a 


P 


Aconf 


z f 


K 


N c 

1 " conr 


z f 


K 


8 


6.0219 


1000 


1.74 
1.77 
1.80 
1.83 
1.86 


0.1534 
0.1537 


1000 


1.79 
1.80 
1.81 
1.82 


0.1530 
0.1534 
0.1537 
0.1540 


10 


6.1628 


1000 


1.73 
1.76 
1.79 
1.82 


0.1521 
0.1522 


1000 


1.78 
1.79 
1.80 
1.81 


0.1520 
0.1521 
0.1522 
0.1523 


12 


6.2885 


500 


1.71 
1.74 
1.77 
1.80 


0.15050 
0.15100 


300 


1.70 
1.73 
1.77 
1.80 


0.15025 
0.15050 
0.15100 
0.15125 


16 


6.4956 


300 


1.64 
1.67 
1.70 
1.73 
1.76 


0.1489 
0.1490 


100 


1.70 
1.71 
1.73 
1.74 


0.1488 
0.1489 
0.1490 
0.1491 


20 


6.6790 


112 


1.66 
1.68 
1.70 
1.72 


0.1473 
0.1474 
0.1475 
0.1476 








24 


6.8187 


100 


1.60 
1.63 
1.66 
1.69 


0.1463 
0.1464 
0.1465 
0.1466 









Table C.l: Values of k and Zf used for the tuning and number of configurations, N con {, used in all calculations at the 
corresponding value of /?. Scale NP (see text). The data of the last three columns have been used only for a separate analysis 
with method (1), which we denote here as method (1*). 
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Guess values for the tuning 
Intermediate scale: g 2 = 2.4484 



L/a 


P 


N c 

1 ' colli 


z f 


K 


N t 

1 " conr 


z f 


K 


8 


7.0197 


1000 


1.51 


0.14445 


1000 


1.35 


0.14440 








1.54 


0.14450 




1.45 


0.14445 








1.57 






1.55 


0.14450 








1.60 






1.65 


0.14455 


12 


7.3551 


500 


1.46 


0.1431 


300 


1.50 


0.1430 








1.49 


0.1432 




1.51 


0.1431 








1.52 






1.52 


0.1432 








1.55 






1.53 


0.1433 


16 


7.6101 


300 


1.44 


0.1421 


100 


1.48 


0.1420 








1.47 


0.1422 




1.49 


0.1421 








1.50 






1.50 


0.1422 








1.53 






1.51 


0.1423 



Table C.2: Same caption as in tab. C.l but at scale I (see text) 





Guess values for the tuning 








Perturbative scale: g 2 = 


= 0.9944 






L/a 





N t 

1 ' colli 




K 


N 

1 " colli 


*/ 


K 


8 


10.3000 








1000 


1.2955 


0.13541 














1.2965 


0.13544 














1.2975 


0.13547 














1.2985 


0.13550 


12 


10.6086 








300 


1.292 


0.13514 














1.294 


0.13517 














1.297 


0.13520 














1.299 


0.13523 


16 


10.8910 


300 


1.23 


0.13484 


100 


1.285 


0.13482 








1.26 


0.13487 




1.286 


0.13484 








1.29 






1.287 


0.13487 








1.32 






1.288 


0.13489 



Table C.3: Same caption as in tab. C.l but at scale P (see text) 
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Guess values for the tuning 




2P scale 






L/a 




N r 
1 " conr 


z f 


K 


16 


12.0000 


100 


1.23 


0.1335 








1.24 


0.1336 








1.25 


0.1337 








1.26 


0.1338 



Table C.4: Same caption as in tab. C.l but at scale 2P (see text). Here no separate tuning was performed for method (1). 





Guess values for the tuning 




PP scale 






L/a 


13 


N e 
1 " conr 


z f 


K 


16 


24.0000 


80 


1.11 


0.1287 








1.12 


0.1288 








1.13 


0.1289 








1.14 


0.1290 



Table C.5: Same caption as in tab. C.l but at scale PP (see text). Here no separate tuning was performed for method (1). 
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C.2 Tuning results 






KorP*) 


Kcr(l) 


K CI (2) 


Kcr (3) 


«cr(4) 


«cr(5) 


K cr (6) 


«cr(7) 


Kor(SP) 








Hadronic scale: L = 


1.436 r (/J ~ 300MeV) 








6.0219 


0.153530 (24) 


0.15353(66) 


0.15354 (66) 


0.15352 (67) 


0.15354(66) 


0.15354 (66) 


0.15354 (66) 


0.15353 (66) 


0.153371 (10) 


6.1628 


0.152134(17) 


0.15213(66) 


0.15214(66) 


0.15213 (67) 


0.15214(66) 


0.15214(66) 


0.15214 (66) 




0.152012 (7) 


6.2885 


0.150815 (22) 


0.15082 (66) 


0.15082 (66) 


0.15082 (66) 


0.15082 (65) 


0.15082 (65) 


0.15082 (65) 


0.15082 (66) 


0.150752 (10) 


6.4956 


0.148945 (25) 


0.14894(34) 


0.14894(33) 


0.14893 (34) 


0.14894(33) 


0.14894(33) 


0.14894 (33) 


0.14894(33) 


0.148876 (13) 


6.6790 




0.14748 (74) 


0.14748 (74) 


0.14748 (74) 


0.14748 (73) 


0.14748 (73) 


0.14748 (73) 






6.8187 




0.14645 (41) 


0.14645 (41) 


0.14645 (42) 


0.14645 (41) 


0.14645 (41) 


0.14645 (41) 


0.14645 (41) 










Intermediate scale: 


9 2 = 2.4484 (fj, ~. 


lGeV) 








7.0197 


0.144501 (13) 


0.14450 (41) 


0.14450 (41) 


0.14450 (41) 


0.14450 (41) 


0.14450 (41) 


0.14450 (41) 


0.14450 (41) 


0.144454 (7) 


7.3551 


0.143113 (12) 


0.14311 (29) 


0.14311 (29) 


0.14311 (29) 


0.14311 (29) 


0.14311 (29) 


0.14311 (29) 


0.14311 (29) 


0.143113 (6) 


7.6101 


0.142112 (13) 


0.14212 (23) 


0.14212 (23) 


0.14212 (23) 


0.14212 (23) 


0.14212 (23) 


0.14212 (23) 


0.14212 (23) 


0.142107 (6) 








Pcrturbative scale: c, 


2 = 0.9944 (m ~ 


30GeV) 








10.3000 


0.1354609 (54) 
















0.135457(5) 


10.6086 


0.1351758 (56) 
















0.135160 (4) 


10.8910 


0.1348440 (61) 


0.134844(93) 


0.134844(93) 


0.134844 (93) 


0.134844 (93) 


0.134844 (93) 


0.134844 (93) 


0.134844(93) 


0.134849 (6) 


2P scale 


12.0000 




0.13363(41) 


0.13363 (41) 


0.13363 (41) 


0.13363 (41) 


0.13363 (41) 


0.13363 (41) 


0.13363 (41) 




PP scale 


24.0000 




0.12877(15) 


0.12877(15) 


0.12877(15) 


0.12877(15) 


0.12877 (15) 


0.12877(15) 


0.12877 (15) 





Table C.6 

Summary table of re cr for all beta values and tuning conditions, (1) to (7) (see sect. |5.2| for a description of all the methods). 
The data of column (1*) correspond to a separate analysis with method (1), using slightly different simulation parameters 
(cf. tab. |C.l| tab. |C.2| and tab. |C.3| . For reference, we also give K cr for the SF I43I6I44| . 



p 


•?(!•) 


*7(D 


*f(2) 


»?(S) 


zf(4) 


-7(5) 


•?(B) 










Hadronic scale: L — 


1.436 t- 








6.0219 


1.8090 (32) 


1.8091 (32) 


1.7946(34) 


1.8434 (37) 


1.7656 (27) 


1.7597 (27) 


1.7835 (40) 


1.7980 (17) 


6.1628 


1.7920 (30) 


1.7923 (29) 


1.7820(31) 


1.8175 (33) 


1.7541 (25) 


1.7497 (25) 


1.7687(38) 




6.2885 


1.7664 (51) 


1.7658 (38) 


1.7573(40) 


1.7869 (46) 


1.7312 (34) 


1.7283 (32) 


1.7408 (56) 


1.7509(22) 


6.4956 


1.7212 (83) 


1.7201 (41) 


1.7132 (44) 


1.7377 (46) 


1.6929(35) 


1.6894(34) 


1.7053 (63) 


1.7076 (21) 


6.6790 




1.6841 (56) 


1.6789(59) 


1.6973 (65) 


1.6582 (52) 


1.6577(52) 


1.6600(90) 




6.8187 




1.6427(56) 


1.6381 (60) 


1.6529 (60) 


1.6253 (51) 


1.6201 (50) 


1.6421 (88) 


1.6366(27) 


Intermediate scale: g 2 — 2.4484 


7.0197 


1.5467 (15) 


1.5404 (16) 


1.5296(17) 


1.5597(18) 


1.5156 (14) 


1.5126 (14) 


1.5229 (21) 


1.5392 (12) 


7.3551 


1.5126 (23) 


1.5139 (18) 


1.5088(19) 


1.5233 (19) 


1.4955 (16) 


1.4945 (15) 


1.4981 (23) 


1.5120 (12) 


7.6101 


1.4942 (37) 


1.4943(20) 


1.4908(21) 


1.5007 (23) 


1.4800(18) 


1.4789 (17) 


1.4827 (28) 


1.4916(13) 


Pcrturbative scale: g 2 — 0.9944 


10.3000 


1.29730 (67) 
















10.6086 


1.2954 (11) 
















10.8910 


1.2858 (15) 


1.28692 (88) 


1.28487 (91) 


1.28984 (99) 


1.27999(83) 


1.27976 (75) 


1.2805 (13) 


1.28619(67) 


2P scale 


12.0000 




1.2493(14) 


1.2481 (14) 


1.2510 (15) 


1.2438 (13) 


1.2443 (12) 


1.2428 (20) 


1.2494 (10) 


PP scale 


24.0000 




1.11268 (57) 


1.11162 (58) 


1.11391 (65) 


1.11013(54) 


1.11003 (48) 


1.11030(87) 


1.11269 (50) 



Table C.7: Summary table of z" for all beta values and tuning conditions, (1) to (7) (see sect. 5.2 for a description of all 
the methods). The data of column (1*) correspond to a separate analysis with method (1), using slightly different simulation 
parameters (cf. tab. IC.ll tab. |C.2|and tab. |C.3||. 
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Azf 


m) for different methods 










Hadronic scale: L 


= 1.436 r 






L/a 


2 


3 


4 


5 


6 


7 


8 


0.0145 (47) 


-0.0343 (49) 


0.0435 (42) 


0.0494 (42) 


0.0256(51) 


0.0111 (36) 


10 


0.0103 (42) 


-0.0252 (44) 


0.0382 (38) 


0.0426 (38) 


0.0236 (48) 




12 


0.0085 (55) 


-0.0211 (60) 


0.0346 (51) 


0.0375 (50) 


0.0250 (68) 


0.0149 (44) 


16 


0.0069 (60) 


-0.0176 (62) 


0.0272 (54) 


0.0307 (53) 


0.0148(75) 


0.0125 (46) 


20 


0.0052 (81) 


-0.0132 (86) 


0.0259 (76) 


0.0264 (76) 


0.024(11) 




24 


0.0046 (82) 


-0.0102 (82) 


0.0174(76) 


0.0226 (75) 


0.001 (10) 


0.0061 (62) 



Table C.8: Differences of z" , Az"(m), determined from different methods. The differences are always z"(l) minus z"{ 
obtained from any other method m = 2 . . . 7. Scale NP. 



Continuum limit of Az^ r (m). 


Hadronic scale: L 


= 1.436 r 




Fit: Azf = b 


+ bif 


m 


bo 


6i 


2 


0.0006(99) 


0.10(12) 


3 


-0.001 (10) 


-0.24(13) 


4 


0.0076(91) 


0.31 (11) 


5 


0.0096(91) 


0.33(11) 


6 


-0.000 (12) 


0.25(15) 


7 


-0.001 (12) 


0.20(18) 



Table C.9: Continuum limit of the data presented in tab. |C.8| We have performed a linear fit in a/L. The point L/a 
is not included in the fit. 
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C.3 Scaling studies 



Tl\' 2 -even correlation functions 
NP scale: L = 1.436 r and = (0, 0, 0) 



L/a 


P 


m 


3p 






M Vo- 


a v 0+ 


si 1 


8 


6.0219 


1 


3.637 


(14) 


0.02609 (25) [82] 


l.9779 (85) 


0.01279 (16) [42] 


1.1857 (53) 






2 


3.643 


(141 


0.02249 (19) [63] 


1.9839 (85) 


0.01090 (12) [32] 


1.1883 (52) 






3 


3.619 


(141 


0.03956 (38) [164] 


1.9568 (87) 


0.01990 (24) [86] 


1.1767 (54) 






4 


3.652 


(13) 


0.01903 (13)[17] 


1.9907(83) 


0.009143 (79) [89] 


1.1915 (51) 






5 


3.653 


(1 31 


01 893 (141[1 41 


1 991 2 ("831 


0091 1 2 ("SO) [801 


1 1 91 8 f 51 1 






6 


3.647 


(13) 


0.02057 (15)[53] 


1.9873 (84) 


0.00991 (10) [27] 


1.1899 (52) 






7 


3.642 


(14) 


0.02322 (20) 


1.9827 (85) 


0.01128 (13) 


1.1878 (52) 


10 


6.1628 


1 


3.605 


(13) 


0.02081 (22) [93] 


1.8788 (84) 


0.01028 (13) [47] 


1.1136 (51) 






2 


3.607 


(131 


0.01738 (18) [85] 


1.8815 (84) 


0.00852 (11) [43] 


1.1143 (51) 






3 


3.598 


(141 


0.03178 (29)[153] 


1.8685 (85) 


0.01592 (17) [78] 


1.1104 (52) 






4 


3.610 


(13) 


0.011182 (99) [312] 


1.8843 (83) 


0.005342 (60) [156] 


1.1144 (50) 






5 


3.610 


(13) 


010616 (88) f260l 


1 8842 (82) 


005053 (53) [1301 


1.1142 (50) 






6 


3.609 


(13) 


0.01386 (14) [80] 


1.8836 (83) 


0.006714 (85) [407] 


1.1147(50) 


12 


6.2885 


1 


3.528 


(191 


0.01809 (24)[129] 


1.808 (12) 


0.00900 (14) [64] 


1.0660 (70) 






2 


3.530 


(191 


0.01510 (22)[118] 


1.810 (11) 


0.00748 (13) [59] 


1.0663 (69) 






3 


3.524 


(191 


0.02724 (30)[210] 


1.801 (12) 


0.01363 (18) [105] 


1.0644 (71) 






4 


3.531 


(181 


0.00836 (13) [54] 


1.811 (11) 


0.004065 (76) [274] 


1.0657 (68) 






5 


3.531 


(18) 


0.00784 (12)[48] 


1.811 (11) 


0.003801 (71)[241] 


1.0655 (68) 






g 


3.531 


(191 


01041 (16M1261 


1.811 (11) 


005104 (94)[638l 


1 0662 (691 






7 


3.530 


(191 


0.01310 (20) 


1.810 (11) 


0.00647(11) 


1.0663 (69) 


16 


6.4956 


1 


3.458 


(22) 


0.01423 (18)[140] 


1.702 (13) 


0.00691 (10) [68] 


0.9921 (77) 






2 


3.458 


(22) 


0.01179 (16)[135] 


1.702 (13) 


0.005714 (90) [658] 


0.9919 (77) 






3 


3.456 


(22) 


0.02169 (22)[209] 


1.698 (13) 


0.01057(13)[101] 


0.9918 (78) 






4 


3.456 


(211 


0.00620 (11)[68] 


1.702 (13) 


0.002974 (59) [331] 


0.9905 (76) 






5 


3.456 




005474 ("991 T5961 


1 702 ( 1 31 


002620 C541 [2901 


9902 C7fi1 






6 


3.457 


(211 


0.00933 (14) [174] 


1.703 (13) 


0.004510 (78) [849] 


0.9915 (77) 






7 


3.457 


(21) 


0.01001 (15) 


1.703 (13) 


0.004842 (82) 


0.9917(77) 


20 


6.6790 


1 


3.377 


(35) 


0.01356 (22) [194] 


1.628 (22) 


0.00650 (13) [93] 


0.948 (13) 






2 


3.377 


'35) 


0.01172 (20)[189] 


1.628 (22) 


0.00561 (12) [91] 


0.948 (13) 






3 


3.377 


(35) 


0.01886 (27) [269] 


1.626(22) 


0.00903 (16) [127] 


0.948 (14) 






4 


3.375 


(34) 


0.00571 (13) [115] 


1.627(22) 


0.002713 (72) [554] 


0.946(13) 






5 


3.375 


(34) 


0.00559 (13) [115] 


1.627(22) 


0.002656 (71) [555] 


0.946 (13) 






6 


3.375 


(34) 


0.00614 (14)[215] 


1.628(22) 


0.002925 (76) [1037] 


0.946(13) 


24 


6.8187 


1 


3.346 


(38) 


0.00805 (14) [150] 


1.643(21) 


0.003932 (81) [728] 


0.950(13) 






2 


3.346 


(37) 


0.00679 (13) [147] 


1.643(21) 


0.003318 (73) [715] 


0.950(13) 






3 


3.345 


(38) 


0.01125 (18)[196] 


1.642(21) 


0.00549 (10) [95] 


0.951 (13) 






4 


3.345 


(37) 


0.003896 (90) [915] 


1.643(21) 


0.001900 (50) [445] 


0.950(13) 



Table C.10: (continuing on the next page, caption below) 
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1 2 

7?- 5 ' -even correlation functions 










NP scale: L = 1.436 r and 6 = (0, 0, 0) 




L/a 




m 




a 11 


a 12 
9 v - 


9v 0+ 


a} 1 






5 


3.344 (37) 


0.002970 (75) [765] 


1.643(21) 


0.001447 (42) [372] 


0.949(13) 






6 


3.346 (38) 


0.00788 (14) [243] 


1.643(21) 


0.003849 (80)[1188] 


0.950(13) 






7 


3.345 (37) 


0.00641 (12) 


1.643(21) 


0.003129 (70) 


0.950(13) 



Table C.IO: Results for the TZ-^-even correlation functions, gp 1 , g^ g and g^ 1 , at scale L = 1.436 ro and for 6 = (0,0,0). 
The data have been obtained via linear interpolations to th e cr itical val ues o f k and quadratic interpolations to the critical 
values of zj for all the tuning methods (1) to (7). See sect. |6.1|and sect. |6.2|for explanations. For the correlation functions 
gp 1 ^ and ffy^ ^ rs * error * s statistical while the second, where available, contains the propagation of the error in z"' 

(sec text). 









7?.g' 2 -even correlation functions 










I scale: g 2 = 2.4484 and 6 = (0, 0, 0) 






L/a 


/3 


m 






<- 


a 12 


9V 


8 


7.0197 


1 


3.2342(83) 


0.010013 (72) [84] 


2.2115(68) 


0.006152 (52) [57] 


1.5371 (48) 






2 


3.2391 (82) 


0.010043 (76) [82] 


2.2170(68) 


0.006203 (54) [62] 


1.5412(48) 






3 


3.2240 (85) 


0.011667 (97) [237] 


2.1987 (69) 


0.007214 (71) [153] 


1.5286 (49) 






4 


3.2446(81) 


0.01110 (10) [13] 


2.2222 (67) 


0.006959 (69) [97] 


1.5456 (47) 






5 


3.2456(81) 


0.01148 (11) [14] 


2.2231 (67) 


0.007222 (74) [102] 


1.5464 (47) 






6 


3.2418(82) 


0.010407 (86) [142] 


2.2197(67) 


0.006468 (59) [104] 


1.5434 (47) 






7 


3.2348 (83) 


0.009982 (72) 


2.2121 (68) 


0.006135(52) 


1.5376 (48) 


12 


7.3551 


1 


3.219(11) 


0.003185 (35) [99] 


2.1342(89) 


0.001976 (25) [64] 


1.4657(61) 






2 


3.220(11) 


0.002968 (32) [58] 


2.1358(89) 


0.001835 (22) [37] 


1.4668 (61) 






3 


3.216(11) 


0.004005 (48) [201] 


2.1305(89) 


0.002514 (35) [131] 


1.4633 (61) 






4 


3.223(10) 


0.003154 (40) [71] 


2.1386(88) 


0.001966 (26) [49] 


1.4691 (61) 






5 


3.224(10) 


0.003212 (42) [77] 


2.1387(88) 


0.002005 (27) [53] 


1.4693 (61) 






6 


3.223(10) 


0.003032 (37) [80] 


2.1382(88) 


0.001884 (24) [55] 


1.4688 (61) 






7 


3.220(11) 


0.003085 (34) 


2.1348(89) 


0.001911 (24) 


1.4662 (61) 


16 


7.6101 


1 


3.203(13) 


0.001786 (30) [155] 


2.094(11) 


0.001123 (21)[100] 


1.4314(76) 






2 


3.204(13) 


0.001526 (25)[123] 


2.095(11) 


0.000954 (18) [80] 


1.4319(75) 






3 


3.202(13) 


0.002454 (40) [256] 


2.092(11) 


0.001559 (28) [166] 


1.4303 (76) 






4 


3.206(13) 


0.001195 (17) [17] 


2.096(11) 


0.000741 (12) [12] 


1.4331 (75) 






5 


3.206(13) 


0.001201 (18)[20] 


2.096(11) 


0.000745 (12) [14] 


1.4332(75) 






6 


3.205(13) 


0.001211 (18)[41] 


2.096(11) 


0.000750 (12) [26] 


1.4329(75) 






7 


3.204(13) 


0.001579(26) 


2.095(11) 


0.000988 (18) 


1.4318(75) 



Table C.ll: Same caption as in tab. C.IO but at scale g J = 2.4484 
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1 2 

7^ 5 ' -even correlation functions 










P scale: g 2 = 0.9944 and 6 = (0, 0, 0) 






L/a 


P 


m 




9l\ 


a 12 


a 12 


91 1 


16 


10.8910 


1 


3.0565 (76) 


0.0006071 (87) [89] 


2.4403 (81) 


0.0004638 (70) [72] 


1.9652(65) 






2 


3.0568(76) 


0.0006251 (90) [148] 


2.4407 (81) 


0.0004782 (72) [121] 


1.9656(65) 






3 


3.0560 (76) 


0.0006260 (90) [164] 


2.4396 (81) 


0.0004787 (73) [129] 


1.9646 (65) 






4 


3.0574(75) 


0.000772 (11)[32] 


2.4414(81) 


0.0005952 (88) [263] 


1.9663 (65) 






5 


3.0574(75) 


0.000782 (11) [30] 


2.4414(81) 


0.0006035 (89) [245] 


1.9663 (65) 






6 


3.0573(75) 


0.000750 (11)[49] 


2.4414(81) 


0.0005775 (86) [395] 


1.9662 (65) 






7 


3.0566(76) 


0.0006105 (88) 


2.4405 (81) 


0.0004665 (70) 


1.9653 (65) 



Table C.12: Same caption as in tab. C.10 but at scale g 



0.9944. 







7?.g' 2 -even correlation functions 












2P scale and = 


'0,0,0) 






L/a 


P 


m 


9l\ 


a 11 
9 P+ 


a 12 


a 12 


9\ l 


16 


12.0000 


1 


3.049(11) 


0.000537 (16) [28] 


2.527(13) 


0.000429 (13) [23] 


2.089(11) 






2 


3.049(11) 


0.000569 (17) [37] 


2.527(13) 


0.000455 (14) [31] 


2.089(11) 






3 


3.048(11) 


0.000508 (15) [20] 


2.526(13) 


0.000404 (12) [17] 


2.088(11) 






4 


3.049(11) 


0.000754 (21) [62] 


2.527(13) 


0.000609 (18) [52] 


2.090(11) 






5 


3.049(11) 


0.000726 (21) [50] 


2.527(13) 


0.000586 (17) [42] 


2.090(11) 






6 


3.049(11) 


0.000813 (23) [112] 


2.527(13) 


0.000658 (19) [94] 


2.090(11) 






7 


3.049(11) 


0.000535(16) 


2.526(13) 


0.000427(13) 


2.089(11) 



Table C.13: Same caption as in tab. C.10 but at scale 2P. 









1 2 

7?- 5 ' -even correlation functions 












PP scale and = 


(0,0,0) 






LI a 


P 


m 


9\\ 




a 12 


o 12 

9y + 


91 1 


16 


24.0000 


1 


3.0201(56) 


0.0002171 (50) [61] 


2.7853 (73) 


0.0001936 (46) [57] 


2.5652(67) 






2 


3.0201 (56) 


0.0002289 (52)[81] 


2.7854 (73) 


0.0002045 (48) [76] 


2.5653(67) 






3 


3.0200 (56) 


0.0002117 (50) [50] 


2.7852 (73) 


0.0001885 (46) [46] 


2.5651 (67) 






4 


3.0202 (56) 


0.0002566 (56)[114] 


2.7855 (73) 


0.0002301 (51) [106] 


2.5655 (67) 






5 


3.0202 (56) 


0.0002589 (56) [99] 


2.7855 (73) 


0.0002322 (51) [92] 


2.5655 (67) 






6 


3.0202 (56) 


0.0002528 (55) [185] 


2.7855 (73) 


0.0002266 (50) [172] 


2.5654 (67) 






7 


3.0201 (56) 


0.0002171 (50) 


2.7853 (73) 


0.0001935 (46) 


2.5652 (67) 



Table C.14: Same caption as in tab. C.10 but at scale PP. 
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1 2 

7£ g ' -even correlation functions 
NP scale: L = 1.436 r and = (0.5,0.5,0.5) 



L /a 


P 




n 


9 P+ 


9v - 


n 12 
9v 0+ 


91 1 


Q 
O 


ft noi n 
O.UZiy 


i 


O O/l 10/ QftA 

Z.y41z (oo ) 


n ri/i Q7fi ( 1 ft A 
U.U4y i (ID) 


1 a 1 on /ft/i \ 
1.41Zy (D4 ) 


n noon7 / 1 o^ 
U.UZZU ( ( 1 Z) 


0.8461 


(39) 






1 


o no 7 (i n\ 

z.yo / (10) 


n n /i n7t> f'oi \ l"7ftl 
U.U4y to (Zl ) |/ 0J 


1.4Uyo (b4) 


n noonft /1 /i\Tq/i1 
U.UzzUd (14)[o4J 


0.8443 


(39) 






2 


2.9421 (98) 


0.04o5o (17)[61J 


1.4150 (do) 


0.02060 (12)[27J 


0.8466 


(39) 






3 


fi^n ( 1 n\ 
z.yzo (10) 


U.UolZo (o0)[14o| 


l.oyzy (do ) 


n no7on / 1 n\ Tftftl 
U.Uz ( 2 y 3 (19) [DOj 


0.8371 


40) 






4 


z.yoiu (yo ) 


U.U4Zyi (10 ) [Z4J 


1.4Z10 (DZy 


n m nnn fi nM"i ol 
u.uiyuu (1UJ[1ZJ 


0.8498 


38) 






c; 
O 


z.yoz^ (yo j 


n 0zi9Rzi (x 1 (\\ ronl 


±.^iZZD (DZ j 


u.uiooy (iu^[iij 


0.8502 


38) 






D 


9 O/l c;G / 0'7^ 

z.y4oy (y f J 


n D/1 /I AO j' 1 A reel 
U.U44Dy (10J [ODJ 


1/11 SO ( flO\ 
1.410Z (DZ j 


U.Uiyfo (llj[Zoj 


f\ n j fin 

0.8480 


39) 






<7 
1 


z.y4uy (yy j 


H D/1 791 /■ 1 fi^ 
U.U4 ( zl (loj 


1/11 /'ftQ\ 
1.41oy (Do ) 


n nonon / 1 Q^ 
u.uzuyu (ioj 


0.8461 


39) 


1U 


A 1 coo 
O.lOZo 


1 


O 001 Q f Q X.\ 

z.yzio (oo ) 


U.UooyZ \i-2) 


1.O40Z (Dl ) 


U.U10 ( Oo ( / o) 


0.7946 


(36) 






1 


2.9199 (94) 


U.UobOO (lo ) |olJ 


1.644:2 (d0) 


o m C7n /i 1 \ foci 

0.01579 (11) [obj 


0.7940 


(36) 






2 


z.yzzo (9o j 


U.U00U0 (lb)[;oJ 


l.o4b/ (dU ) 


n ni a ai a ( c\<\\ Toonl 
0.0144/4 (99)|oz9J 


0.7948 


(36) 






Q 

o 


ni 00 / oftA 

z.yizz (yo j 


U.U40iy (Zo ) [lolj 


1 QQ KA 1 P.\ \ 

l.oo04 (Dl ) 


n mnm ( ~\ a \ r^nl 

u.uiyyi (14) [oyj 


0.7910 


(37) 






/I 
4 


z.yzo ( (yi ) 


U.UZ (00 (1U) [oZJ 


1 q c;nft 1 t;o\ 
l.v5oUu (oy ) 


n monnft f' r 7i\ri/iil 
U.U1ZUUD ( / 1) [141J 


0.7958 


35) 









9 G971 fon^ 

z.yz < 1 (yu ) 


u.uz ( uou (yo ) [z ( oj 


i.oouo (oy ) 


n ni 1 7ft 1 f ftQ \ T"! 9 cr 1 
U.U11/D1 (DoJ[lZ0J 


0.7958 


35) 






D 


z.yzio (yz ) 


U.UOUUZ ( lo j [ 1 oj 


l.o^tyi (oy ) 


u.uiouyy (ooj[ozoj 


0.7955 


36) 


1Z 


ft oost; 
O.ZooO 


1 


Z.oO / (lo ) 


U.UZ044 (lo) 


l.oU4o (oU ) 


n ni oc;f;i /'qq\ 
U.U1Z001 (oo ) 


0.7656 


(48) 






1 


O Q C7 / 1 A ~\ 
2.O0 I (14) 


U.UzoZU (ZU)[11UJ 


1 on c; o / Qf\\ 

l.zy5o (oU ) 


n m oqo /1 o\ f/iol 
U.Olzoy (lz) [4oJ 


0.7602 


(47) 






2 


Z.859 (lo) 


O.U25o9 (lo ) [1U1J 


1.29d9 (79) 


0.01127 (11) [45J 


0.7606 


(47) 






Q 

o 


Z.oOZ (14 ) 


u.Uoo / y(zo)[i / / j 


i.zoy4 (oi ) 


n ni K77 / 1 a \ \i(W 

u.uio ( ( (14) [/ yj 


0.7585 


48) 






4 


9 Qft1 fl Q\ 
Z.OOl (lo ) 


n m 009 / 1 9^ r^i 1 
u.uiyyz (izj[oij 


l.zyy ( y(o) 


n nOQftO/l f'7^M'097l 
U.UUoDy4 {I 0) [ZZ l \ 


0.7610 


46) 









Z.oOZ (lo ) 


u.uiy4o (izj [40J 


1 9nG2 fTQ\ 

i.zyyo \(o) 


n nnQ/iQfi ^70\ronil 

U.UUo40o ( / Z ) [ZU1J 


0.7609 


46) 






ft 



Z.oOl (lo) 


U.UZ1 1 U (14) [1 11J 


1 onm Itc\\ 

i.zyyi ( / y) 


n nno/im /'Qft\r/inftl 

u.uuy4yi (oo)[4yDj 


0.7610 


(47) 






7 
/ 


Z.OOU (lo ) 


U.UZ4UU (1 ( ) 


1 0070 I TQ\ 

i.zy ( y ( i y j 


n m n^i n fcm\ 
u.uiuoiy (yo ) 


0.7609 


(47) 


1 ft 

lo 


0.4y00 


1 


Z.ol / (lo ) 


u.uiyo i (11) 


1 OQ 1 7 /nQ\ 

i.zoi / (yo ) 


U.UUo4oU (Do ) 


0.7138 


55) 






1 


z.oi4 (lo ) 


n nioofi/iftAriiftl 

u.uiyzo (Id) [i loj 


1 0Q111 /no\ 

l.zoui (yz ) 


n nnQOQQ ( Qn\r /in7l 
u.uuoZoo (oy)|4y(j 


0.7129 


(55) 






2 


2.815 (15) 


O.UlYZo (14)[llzJ 


"i oun /oo\ 

l.zolO (9z) 


0.007oo4 (ol)[4o0J 


0.7129 


55) 






q 
O 


Z.olZ (lo ) 


n no^/in ri 7/il 
U.UZ04U (iyj [1 (4J 


1 oo^a I qi\ 
i.zzdo (yo ) 


u.uiuyu (11 ) [ ( oj 


0.7122 


55) 






4 


Z.olO (lo ) 


u.uizooo (yyj[oyuj 


1 9Q9n f Q r )\ 

i.zozu (yz ) 


n nnt^QQfi ^nATo^ol 
u.uuoooo (oyj[zozj 


0.7125 


55) 









Z.olo (lo ) 


u.uiiy4u (yoj[oioj 


1 9Q 1 n f Ci r )\ 

i.zoiy (yz ) 


n nn^n7i ^ftAToool 

U.UUOU i 1 (00 ) [ZZZ\ 


0.7124 


55) 






ft 


Z.olO (lo ) 


U.Uloiy (1ZJ [14DJ 


1 9Q 1 ft f Cil\ 

i.zoid (yz ) 


n nnft/i7fi f79M"ftQil 

U.UUD4 ( o ( / Z ) [DolJ 


0.7129 


55) 






I 


Z.olO (lo ) 


n m K7c /i q\ 
U.U10 (0 (lo) 


1 oq 1 £ ^no\ 

l.zoio (yz j 


n nnft70Q C7v^\ 
U.UUD (26 ( / 0) 


0.7129 


55) 


20 


6.6790 


1 


2.709 (24) 


0.01578 (18) [IdOJ 


1.16b (15) 


U.0U672 (11) [d7J 


0.6748 


(92) 






2 


2.709 (24) 


0.01427 (16) [156] 


1.166(15) 


0.00608 (10) [65] 


0.6747 


(92) 






3 


2.708 (24) 


0.02009 (21) [222] 


1.164(15) 


0.00856 (14) [94] 


0.6746 


(92) 






4 


2.709 (23) 


0.00934 (11) [97] 


1.167(15) 


0.003954 (70) [413] 


0.6741 


(92) 






5 


2.709 (23) 


0.00924 (11) [97] 


1.167(15) 


0.003911 (69)[413] 


0.6741 


(92) 






6 


2.709 (23) 


0.00970 (12) [180] 


1.167(15) 


0.004110 (73) [770] 


0.6742 


(92) 



Table C.15: (continuing on the next page, caption below) 
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1 2 

7?- 5 ' -even correlation functions 










NP scale: L = 1.436 ro and 9 = (0.5,0.5,0.5) 




L/a 


R 
P 


m 


9 P _ 




n 12 


a 12 
9v 0+ 


y\ 


OA 
Z4 


u.olo f 


1 


Z.D / O \Z i ) 


U.UUy /b (loj [IzzJ 


l.lOD ^IbJ 


n nri/ii c\i ('7^\\z^~\ ol 
U.UU411U (/Oj|OlzJ 


n ^90 (Q1\ 






2 


2.D/5 (27) 


0.008 to (12) [120J 


1.13b (lb) 


U.UUobfo ( i Oo)[oUUJ 


^90 (CM\ 






*? 
o 


2.675 (27) 


0.01234 (16) [159] 


1.135(16) 


0.005194 (92) [678] 


O R^OR (QA\ 
u.oozo i^y^t J 






A 


2.675 (27) 


0.006389 (88) [756] 


1.136(16) 


0.002683 (50)[317] 


u.dozd i^yo ) 






5 


2.675 (27) 


0.005632 (76) [632] 


1.136(16) 


0.002363 (44) [264] 


0.6524 (93) 






6 


2.675 (27) 


0.00963 (13) [198] 


1.136(16) 


0.004049 (74) [839] 


0.6529 (93) 






7 


2.675 (27) 


0.00843(12) 


1.136(16) 


0.003546 (66) 


0.6529 (93) 



Table C.15: Same caption as in tab. |C,10|but at 9 = (0.5. 0.5,0 .5). The values denoted with (1*) have been obtained by 



direct simulations at the critical values of 
and ~ 12 
text) 



and Zf. See sect. 



5.3 



9v (j+ the first error is statistical while the second, where available, contains the propagation of the error in z" (see 



for explanation of (1*). For the correlation functions gip 

cr 
/ 











7?.g' 2 -even correlation functions 












I scale: g 2 = 2.4484 and 9 = (0.5,0.5,0.5) 




L/a 





m 




yp + 


a v _ 


g\? 


9} 1 


8 


7.0197 


1* 


Z.014Z 


f AQ\ 
4o ) 


003341 f79l 


1 4348 (43) 


01 7953 (661 


1.0006 (30) 






1 


2.5143 


55) 


003306 f 791 11 191 


1 4358 (441 


01 7784 (6511761 


1.0013 (31) 






2 


2.5186 


(54) 


0.003279 (86) [86] 


1.4399 (44) 


0.017680 (66) [66] 


1.0042 (31) 






3 


2.5056 


(56) 


0.003480 (85) [247] 


1.4265(45) 


0.018668 (71) [139] 


0.9953 (31) 






4 


2.5235 


(53) 


0.00332 (10) [90] 


1.4442 (43) 


0.017964 (72) [72] 


1.0075 (31) 






5 


2.5245 


(53) 


0.00334 (11) [93] 


1.4450 (43) 


0.018086 (74) [74] 


1.0081 (30) 






6 


2.5210 


(54) 


0.003288 (93) [93] 


1.4421 (44) 


0.017757 (68) [71] 


1.0058 (31) 






7 


2.5148 


(55) 


0.003300 (79) 


1.4363 (44) 


0.017759 (65) 


1.0016(31) 


12 


7.3551 


1* 


2.5015 


(62) 


0.014097(43) 


1.3856 (53) 


0.007547 (35) 


0.9509 (37) 






1 


2.5003 


(65) 


0.014158 (43)[109] 


1.3847(52) 


0.007568 (36) [60] 


0.9504 (36) 






2 


2.5016 


(65) 


0.013931 (42) [79] 


1.3859 (52) 


0.007449 (35) [45] 


0.9512(36) 






3 


2.4976 


(66) 


0.014895 (47) [190] 


1.3819(53) 


0.007964 (40) [103] 


0.9486 (36) 






4 


2.5044 


(65) 


0.013914 (50) [50] 


1.3883 (52) 


0.007454 (36) [40] 


0.9530 (36) 






5 


2.5046 


(65) 


0.013947 (51)[51] 


1.3884(52) 


0.007473 (36) [40] 


0.9531 (36) 






6 


2.5039 


(65) 


0.013852 (48) [51] 


1.3879(52) 


0.007417 (35) [41] 


0.9527(36) 






7 


2.5008 


(65) 


0.014059 (42) 


1.3852(52) 


0.007516(36) 


0.9507(36) 


16 


7.6101 


1* 


2.4864 


(82) 


0.008178(31) 


1.3483 (67) 


0.004364 (25) 


0.9199 (46) 






1 


2.4868 


(84) 


0.008141 (33) [140] 


1.3484(67) 


0.004335 (27) [75] 


0.9200 (46) 






2 


2.4875 


(84) 


0.007918 (31)[115] 


1.3490 (67) 


0.004217 (26) [61] 


0.9204 (46) 






3 


2.4855 


(85) 


0.008695 (38) [222] 


1.3469 (67) 


0.004631 (30) [120] 


0.9191 (46) 






4 


2.4893 


(84) 


0.007593 (29) [32] 


1.3505 (67) 


0.004047 (23) [23] 


0.9214(46) 






5 


2.4895 


(84) 


0.007591 (29) [29] 


1.3506 (67) 


0.004047 (23) [23] 


0.9215(46) 






6 


2.4889 


(84) 


0.007623 (29) [59] 


1.3502 (67) 


0.004062 (23) [34] 


0.9212(46) 






7 


2.4874 


(84) 


0.007964(31) 


1.3489 (67) 


0.004241 (26) 


0.9203 (46) 



Table C.16: Same caption as in tab. C.15 but at scale g 
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1 2 

7?- 5 ' -even correlation functions 










P scale: g 2 = 0.9944 and = (0.5,0.5,0.5) 




L/a 


P 


m 




9l\ 




a 12 


si 1 


16 


10.8910 


1* 


2.2924 (43) 


0.006800(15) 


1.4342 (44) 


0.004191 (14) 


1.1555 (35) 






1 


2.2921 (44) 


0.006796 (15) [16] 


1.4340 (44) 


0.004181 (14) [14] 


1.1553 (35) 






2 


2.2924 (44) 


0.006796 (16) [16] 


1.4343 (44) 


0.004182 (14) [14] 


1.1556 (35) 






3 


2.2916(45) 


0.006830 (15) [25] 


1.4334 (44) 


0.004202 (14) [18] 


1.1549 (35) 






4 


2.2931 (44) 


0.006872 (17) [21] 


1.4350 (44) 


0.004230 (15) [18] 


1.1562 (35) 






5 


2.2931 (44) 


0.006878 (17) [21] 


1.4350 (44) 


0.004234 (15) [18] 


1.1563 (35) 






6 


2.2930 (44) 


0.006859 (17) [29] 


1.4349 (44) 


0.004222 (15) [22] 


1.1562 (35) 






7 


2.2922 (44) 


0.006794(15) 


1.4341 (44) 


0.004180(14) 


1.1554 (35) 



Table C.17: Same caption as in tab. C.15 but at scale g 



0.9944. 







7£g' 2 -even correlation functions 












2P scale and 9 = (0.5,0.5,0.5) 






L/a 


P 


m 




a 11 


a 12 


a 12 
9v 0+ 


'A 1 


16 


12.0000 


1 


2.2686 (65) 


0.006722 (22) [22] 


1.4486 (62) 


0.004247 (19) [20] 


1.1988(52) 






2 


2.2688 (65) 


0.006737 (22) [25] 


1.4488 (62) 


0.004257 (19) [23] 


1.1990(52) 






3 


2.2683 (66) 


0.006712 (21) [21] 


1.4483 (62) 


0.004240 (19) [19] 


1.1985(52) 






4 


2.2694 (65) 


0.006845 (25) [40] 


1.4494 (62) 


0.004327 (21) [31] 


1.1996(51) 






5 


2.2693 (65) 


0.006828 (25) [33] 


1.4493 (62) 


0.004316 (20) [27] 


1.1995(51) 






6 


2.2695 (65) 


0.006882 (26) [71] 


1.4495 (62) 


0.004351 (21) [50] 


1.1997(51) 






7 


2.2686 (65) 


0.006721 (22) 


1.4486 (62) 


0.004246(19) 


1.1988(52) 



Table C.18: Same caption as in tab. C.10 but at scale 2P and B = (0.5,0.5,0.5) 







7?.g' 2 -even correlation functions 












PP scale and = (0.5,0.5,0.5) 






L/a 


P 


m 




a 11 
9 F+ 


a 12 


a 12 


9Y 


16 


24.0000 


1 


2.1966(29) 


0.006459 (11) [11] 


1.5009 (31) 


0.004391 (10) [10] 


1.3851 (28) 






2 


2.1967(29) 


0.006459 (11) [11] 


1.5011 (31) 


0.004392 (10) [11] 


1.3853 (28) 






3 


2.1964(29) 


0.006464 (10) [11] 


1.5007(31) 


0.004395 (11) [11] 


1.3850(28) 






4 


2.1968(29) 


0.006468 (11) [11] 


1.5013(31) 


0.004398 (10) [11] 


1.3855(28) 






5 


2.1969(29) 


0.006469 (11) [11] 


1.5013(31) 


0.004399 (10) [11] 


1.3855(28) 






6 


2.1968(29) 


0.006467 (11) [13] 


1.5012(31) 


0.004397 (10) [12] 


1.3854(28) 






7 


2.1966(29) 


0.006459(11) 


1.5009(31) 


0.004391 (10) 


1.3851 (28) 



Table C.19: Same caption as in tab. C.10 but at scale PP and 6 = (0.5,0.5,0.5) 



44 



7£c' 2 -odd correlation functions 


NP scale: L = 1.436 7*0 and = 


= (0, 0, 0) 


T In 


a 
P 


m 




77 12 


8 


6.0219 


1 




-U.171 (1U ) 






2 


U.U434 (97) 


n i 1 a /in 1 ! 
-U.114 (1UJ 






3 


-U.lUol (95) 


n one ^ 1 rA 

-U.oUo (1U) 






A 


U.loUO (yes ) 











u. i^toz ^yo j 


o 09^ i "\ 

U.UZO ^ 1 1 ) 






D 


O 07R7 f 07^ 

u.u id / t^y / ^ 


n 070 f 

-U.U IU ^ 1 1 J 






( 


U.Uoo2 {J 1 ) 


o 1 97 fi 

-U. 1Z I ^ lUj 


10 


6.1628 


1 




n 1 A A O. f OQ\ 

-U.i44y (yo ) 






2 


U.U29o (84) 


n 1 nro /no N 

-U.lUoo (9o) 






3 


-U.Uf24 (o4) 


n o /i nr. /no\ 
-U.24U5 (92 ) 






A 


U.iuyo (00 ) 











U.1224 (00 ) 


01 R7 fCt^ 

u.uio / ^yo ) 









U.UO (0 (00 J 


-U.U004 (y4 j 


12 


6.2885 


1 




-U.1Z4 (12) 






2 


U.U2o (11) 


n noQ / 1 o^ 
-U.U9o (12) 






3 


-U.Uoo (1U) 


n i on ( 1 oA 
-U.199 (12) 






A 
4 


U.Uy4 (11; 








5 


U.1U2 (11) 


n m fi /1 oA 
U.U1U (12J 






6 


U.Uoo (11) 


n nQ/i / 1 oA 
-U.Uo4 (12J 






7 


n A/11 / 1 1 \ 
U.U41 (11) 


-U.U ( 1 (12 ) 


16 


o.49oo 


1 




-U.U95 (12 ) 






2 


U.Ulo (11) 


-U.U ( 1 (12 ) 






3 


-U.U4 J (11) 


-U.loo (12) 






/i 
*± 


079 h 1^ 








c 



U.UoZ ^ LI) 


nni9fi9'i 

U.U-LZ ^ LLj 






D 


0^0, hi^ 
u.uoy ^ 1 1) 


~U.U4:0 \ i-Z,J 






/ 


n oqq /" 1 1 A 
U.Uoo (11) 


-U.UOl \ i-^J 


20 


0.079U 


1 




-U.Uo4 (1 / ) 






2 


n m q f 1 A\ 
U.Ulo (14) 


-U.U67 (17) 






q 



-U.Uoo (14 ) 


-U.12o (loj 








U.Uoo (14 J 








5 


U.Uoo (14) 


009 ^17^ 
U.UUZ \i- 1 j 






D 


U.UOl 1, 1^ ) 


-U.UUo (1 / ) 


24 


6.8187 


1 




O O^fi /'I 7^ 

-U.Uoo (1 / ) 






2 


0.012(15) 


-U.U4o (1 / ) 






3 


-0.027(15) 


-0.092 (17) 






4 


0.045(15) 








5 


0.059(15) 


0.017(17) 






6 


0.001 (15) 


-0.056 (17) 



Table C.20: (continuing on the next page, caption below) 
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7?.g' 2 -odd correlation functions 
NP scale: L = 1.436 r and 6 = (0, 0, 0) 


L/a 


P 


m 
7 


0.016(15) 


-0.038(17) 







Table C.20: Results for the 7?.g' 2 -odd correlation functions, gj^ g and gXf h , at scale L = 1.436 ro and for = (0,0,0). 
The data have been obtained via linear interpolations to the critical values of re and zj for all the tuning methods (1) to 
(7). The lines represent the cases where the corresponding quantity has been used as tuning condition. See sect. |6.3| for 
explanations. 



46 



-odd correlation functions 
I scale: g 2 = 2.4484 and = (0, 0, 0) 



L/a 



a 11 

9 A _ 



77 1 - 



7.0197 



0.0406 (61) 
-0.0730 (60) 
0.0935 (62) 
0.1048 (63) 
0.0659 (62) 
0.0044 (61) 



-0.1134(64) 
-0.0640 (65) 
-0.2015(62) 

0.0137(66) 
-0.0334(65) 
-0.1079 (64) 



12 



7.3551 



0.0192(67) 
-0.0355 (67) 
0.0695 (68) 
0.0732 (68) 
0.0597 (68) 
0.0072 (67) 



-0.0845 (72) 
-0.0610 (72) 
-0.1277(71) 



0.0048 (73) 
-0.0118(72) 
-0.0757(72) 



16 



7.6101 



0.0131 (75) 
-0.0236 (74) 
0.0531 (75) 
0.0572(75) 
0.0431 (75) 
0.0101 (75) 



-0.0645 (80) 
-0.0487(80) 
-0.0933 (80) 



0.0048 (81) 
-0.0123(80) 
-0.0523 (80) 



Table C.21: Same caption as in tab. C.20 but at scale g 2 = 2.4484 



7£r' 2 -odd correlation functions 

5 


P scale: g 2 = 0.9944 and = 


(0,0,0) 


L/a 


P 


m 


9 A _ 


a 12 


16 


10.8910 


1 




-0.0363 (43) 






2 


0.0096 (41) 


-0.0255 (43) 






3 


-0.0136 (41) 


-0.0516 (43) 






4 


0.0323 (41) 








5 


0.0334 (41) 


0.0012 (44) 






6 


0.0299 (41) 


-0.0027 (44) 






7 


0.0034 (41) 


-0.0325 (43) 



Table C.22: Same caption as in tab. C.20 but at scale Y = 0.9944 



47 



1 2 

7?- 5 ' -odd correlation functions 


2P E 


calc and = (0, 0, 0) 


L/a 


P 


m 


9a _ 




16 


12.0000 


1 




-0.0292 (70) 






2 


0.0059 (67) 


-0.0228 (70) 






3 


-0.0081 (67) 


-0.0382 (70) 






4 


0.0266 (67) 








5 


0.0242 (67) 


-0.0026 (70) 






6 


0.0315(67) 


0.0054(70) 






7 


-0.0004 (67) 


-0.0297(70) 



Table C.23: Same caption as in tab. C.20 but at scale 2P. 



PP scale and 6 = (0, 0, 0) 



L/a 



P 



g A _ 



16 



24.0000 



0.0057(31) 
-0.0066 (31) 
0.0137(31) 
0.0143 (31) 
0.0128 (31) 
-0.0001 (31) 



-0.0144(31) 
-0.0084 (31) 
-0.0213(31) 



0.0005 (31) 
-0.0010(31) 
-0.0144 (31) 



Table C.24: Same caption as in tab. C.20 but at scale PP. 



72.ji' 2 -odd correlation functions 
NP scale: L = 1.436 r and 6 = (0.5, 0.5, 0.5) 



L/a 


P 


m 




o 12 


8 


6.0219 


1* 


-0.0270 (59) 








1 


-0.0307(71) 


-0.1465(75) 






2 




-0.1034(76) 






3 


-0.1031 (70) 


-0.2482 (73) 






4 


0.0612 (73) 


-0.0174(78) 






5 
6 


0.0736 (73) 
0.0234 (72) 


-0.0705 (77) 






7 


-0.0072 (71) 


-0.1135(76) 


10 


6.1628 


1* 


-0.0173 (50) 








1 


-0.0210(64) 


-0.1239(70) 






2 




-0.0939 (70) 






3 


-0.0724 (64) 


-0.1971 (69) 






4 


0.0570 (65) 


-0.0129(71) 



Table C.25: (continuing on the next page, caption below) 
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7?.g' 2 -odd correlation functions 


NP scale: L = 


= 1.436 r and 9 = (0.5,0.5,0.5) 


L/a 


a 
P 















U.Uooy (oo ) 








D 


n novo /'fi/i \ 
U.Uz ( Z (04 J 


n nc; rrq /''7n\ 
-U.UoOo ( i V) 


12 


6.2885 


1 


n no i k r: o \ 

-U.UzlD (Do] 








1 


-U.Uluo \( () 


-U.lUoo {ot) 






2 




-U.UoUU (57; 






Q 

o 


-u.uo to [t ( ) 


-U.lol ( (oo ) 






A 
4 


U.UoUo {to) 


-U.UU ( y (oy j 






r: 



U.UD04 ( / o) 








D 


U.Uozl {to) 


n HQ/l /l /eo\ 
-U.Uo44 (oo ) 






I 


U.Ulzo ( / o) 


-U.Uozo (oo ) 


16 


6.4956 


1 


-U.UUyz (5U) 








1 


n m qo f'o yi "\ 
-U.UloZ (o4j 


n no q T /ni \ 

-U.Uoo t (yi j 






2 




-U.Uo4o (91) 






Q 
O 


-U.U4 / 1 (o4 ) 


-u.ioio (yi j 








u.uoyi {o*± ) 


-u.uuyo (yz j 









U.ULOO (O^t j 








D 


n m £0 (qa \ 
U.UlDz (o4 J 


-U.U400 (yz j 






7 


c\ n 1 no f q a\ 
U.UlUo (o4j 


-u.U4yo (yz j 


20 


6.679U 


1 


-U.UU9 (1UJ 


-U.UbY {16) 






2 




-U.U54 (loj 






O 


-U.Uoo (1U J 


n t nn Q^ 
-U.1UU \i~o) 






A 

4 


U.Uo / (11 ) 


n nm ( i Q^ 
-U.UU1 (^"Jj 









n noo. f 1 1 
U.Uoo (11,) 








D 


U.Uo4 (11; 


n nnfi / i q\ 
-U.UUd \ i-o) 


24 


6.8187 


1 


n nno /i 1 \ 
-U.UUo (HJ 


r\ nc:7 / i o\ 
-U.U5 ( (loj 






2 




-U.U4o (lo; 






3 


-0.027(11) 


-0.083(13) 






4 


0.023(11) 


-0.013(13) 






5 


0.032(11) 








6 


-0.007(11) 


-0.056(13) 






7 


0.003(11) 


-0.042(13) 



Table C.25: Same caption as in tab. |C.20| but at = (0^5, 0.5, 0.5). The values denoted with (1*) have been obtained by 
direct simulations at the critical values of k and Zf. See sect. |5.3| for explanation of (1*). 
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-odd correlation functions 
I scale: g 2 = 2.4484 and = (0.5,0.5,0.5) 



L/a 



a 11 

9 A - 



77 1 - 



7.0197 



12 



7.3551 



-0.0378 (31) 
-0.0263 (41) 

-0.0730 (40) 
0.0338 (42) 
0.0411 (42) 
0.0161 (41) 
-0.0234 (41) 
-0.0083 (31) 
-0.0125 (47) 



-0.0355 (47) 
0.0324 (47) 
0.0348 (47) 
0.0260 (47) 

-0.0079 (47) 



-0.0898 (43) 
-0.0550 (43) 
-0.1521 (41) 
-0.0098 (44) 

-0.0334 (44) 
-0.0860 (43) 

-0.0640 (50) 
-0.0472 (50) 
-0.0951 (49) 
-0.0032 (50) 



-0.0118(50) 
-0.0577(50) 



16 



7.6101 



-0.0137(30) 
-0.0084 (50) 

-0.0236 (50) 
0.0257(50) 
0.0283 (50) 
0.0193(50) 

-0.0019(50) 



-0.0500 (54) 
-0.0386 (54) 
-0.0708 (54) 
-0.0035 (54) 

-0.0123(54) 
-0.0412(54) 



Table C.26: Same caption as in tab. C.25 but at scale g z = 2.4484 
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Tt^-odA correlation functions 


P scale: g 2 = 


= 0.9944 and = (0.5, 0.5, 0.5) 


L/a 


P 


m 




a 12 


8 


10.3000 


1* 


-0.0146(14) 




12 


10.6086 


1* 


-0.0063(14) 




16 


10.8910 


1* 


-0.0034(15) 








1 


-0.0056(25) 


-0.0256 (27) 






2 




-0.0183 (27) 






3 


-0.0136(25) 


-0.0360 (27) 






4 


0.0134(25) 


-0.0008 (27) 






5 


0.0140(25) 








6 


0.0120(25) 


-0.0027(27) 






7 


-0.0036 (25) 


-0.0230 (27) 



Table C.27: Same caption as in tab. C.25 but at scale g 2 = 0.9944 



7£g' 2 -odd correlation functions 


2P scale and 9 = (0.5,0.5,0.5) 


L/a 


P 


m 






16 


12.0000 


1 


-0.0034 (40) 


-0.0178 (42) 






2 




-0.0135(42) 






3 


-0.0081 (40) 


-0.0239 (42) 






4 


0.0118(40) 


0.0019 (42) 






5 


0.0104 (40) 








6 


0.0146 (40) 


0.0055 (42) 






7 


-0.0037 (40) 


-0.0181 (42) 



Table C.28: Same caption as in tab. C.20 but at scale 2P and 6 = (0.5,0.5,0.5) 



Tll^-odd correlation functions 

o 


PP scale and = (0.5,0.5,0.5) 


L/a 


P 


m 


9 A _ 


a 12 


16 


24.0000 


1 


-0.0031 (17) 


-0.0097(17) 






2 




-0.0058 (17) 






3 


-0.0066 (17) 


-0.0141 (17) 






4 


0.0043 (17) 


-0.0004 (17) 






5 


0.0046 (17) 








6 


0.0038 (17) 


-0.0010(17) 






7 


-0.0031 (17) 


-0.0097(17) 



Table C.29: Same caption as in tab. C.20 but at scale PP and 6 = (0.5,0.5,0.5) 
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Continuum limit of . NP scale 




e = (o,o,o) 


= (0.5,0.5,0.5) 


m 


bo 


6i 


bo 




1 






0.001 (14) 


-0.22(17) 


2 


-0.001 (18) 


0.31 (23) 






3 


0.003(18) 


-0.75(23) 


0.004(14) 


-0.76(17) 


4 


0.008(18) 


1.02(23) 


0.007(14) 


0.51 (17) 


5 


0.013(18) 


1.09(23) 


0.010(14) 


0.56(17) 


6 


-0.006 (18) 


0.79 (23) 


-0.003(14) 


0.34(17) 


7 


-0.005(29) 


0.57(42) 


-0.005(21) 


0.22(31) 



Table C.30: Continuum limit of for all definitions of z" from (1) to (7). Scale NP. Results are shown for = (0, 0, 0) 

and 8 = (0.5,0.5,0.5), as obtained from linear fits of the data in tab. |C.2"o| and tab. |C25] respectively. The fits are linear 
in a/L and the point L/a = 8 is excluded in all cases. The lines represent the cases where the corresponding quantity has 
been used as tuning condition. See sect. 16.31 for more explanations. 



Continuum limit of g]^ . NP scale 



e = (o,o,o) 

bo &i 



9 = (0.5,0.5,0.5) 
bo &i 



-0.008(21) 
-0.011 (21) 
-0.002 (20) 



-1.37(26) 
-0.96 (26) 
-2.39(26) 



-0.012(16) 
-0.013(16) 
-0.008 (16) 
-0.004(16) 



0.005(21) 
-0.018(21) 
-0.002 (33) 



0.10(26) 
-0.32(26) 
-0.81 (48) 



-0.020(16) 
-0.019(25) 



-1.12(20) 
-0.81 (20) 
-1.89(19) 
-0.07 (20) 



-0.31 (20) 
-0.51 (36) 



Table C.31: Same caption as in tab. C.30 but for 
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